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Abstract 

We establish the behavior of the energy of minimizers of non-local 
Ginzburg-Landau energies with Coulomb repulsion in two space dimen- 
sions near the onset of multi-droplet patterns. Under suitable scaling of 
the background charge density with vanishing surface tension the non- 
local Ginzburg-Landau energy becomes asymptotically equivalent to a 
sharp interface energy with screened Coulomb interaction. Near the on- 
set the minimizers of the sharp interface energy consist of nearly identical 
circular droplets of small size separated by large distances. In the limit 
the droplets become uniformly distributed throughout the domain. The 
precise asymptotic limits of the bifurcation threshold, the minimal energy, 
the droplet radii, and the droplet density are obtained. 

1 Introduction 

Spatial patterns are often a result of the competition between thermodynamic 
forces operating on different length scales. When short-range attractive interac- 
tions are present in a system, phase separation phenomena can be observed, 
resulting in aggregation of particles or formation of droplets of new phase, 
which evolve into macroscopically large domains via coarsening or nucleation 
and growth (see e.g. [1]). This process, however, can be frustrated in the pres- 
ence of long-range repulsive forces. As the droplets grow, the contribution of 
the long-range interaction may overcome the short-range forces, whereby sup- 
pressing further growth. This mechanism was identified in many energy-driven 
pattern forming systems of different physical nature, such as various types of 
ferromagnetic systems, type-I superconductors, Langmuir layers, multiple poly- 
mer systems, etc., just to name a few [2-11]. Remarkably, these systems often 
exhibit very similar pattern formation behaviors [10,12]. 
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One important class of systems with competing interactions are systems in 
which the long-range repulsive forces are of Coulomb type (for an overview, 
see [13, 14] and references therein). The nature of the Coulombic forces may be 
very different from system to system. For example, these forces may arise when 
particles undergoing phase separation carry net electric charge [15-18], or they 
may be a consequence of entropic effects associated with chain conformations 
in polymer systems [19-23]. Coulomb interactions may also arise indirectly 
as a result of diffusion- mediated processes [4,24,25]. All this makes systems 
with repulsive Coulombic interactions a ubiquitous example of pattern forming 
systems. 

Studies of systems with competing short-range attractive interactions and 
long-range repulsive Coulomb interactions go back to the work of Ohta and 
Kawasaki, who proposed a non-local extension of the Ginzburg-Landau energy 
in the context of diblock copolymer systems [19]. Even though its validity for 
diblock copolymer systems may be questioned [21,26-28], the Ohta-Kawasaki 
model is applicable to a great number of physical problems of different ori- 
gin [14]. On the other hand, mathematically Ohta-Kawasaki model presents a 
paradigm of energy-driven pattern forming systems which has been receiving a 
growing degree of attention [9,29-36]. 

The Ohta-Kawasaki energy is a functional of the form [13, 14, 19, 24, 37]: 

£[u] = J (^-\Vu\ 2 + W(u) ) j dx 

+ o / I (u(x) - u)G (x,y)(u(y) - u)dxdy. (1.1) 
z Jn Jn 

Here, u : ft — > R is a scalar quantity denoting the "order parameter" in a 
bounded domain il C M. d . Different terms of the energy are as follows: the 
first term penalizes spatial variations of u on the scales shorter than e, the 
second term, in which W is a symmetric double-well potential drives local phase 
separation towards the minima of IF at u = ±1, and the last term is the long- 
range interaction, whose Coulombic nature comes from the fact that the kernel 
Go solves the Neumann problem for 

-AG (x,y) = S(x-y)-^-, J G (x, y)dx = 0, (1.2) 

where A is the Laplacian in x and 5{x) is the Dirac delta-function. The pa- 
rameter u denotes the prescribed uniform background charge, and the overall 
"charge neutrality" is ensured via the constraint 

— — ■ / udx — u. (1-3) 



It is important to note that the kernel Go solves ( 1.2 1 in the space of the same 
dimensionality as the order parameter u (not to be confused with the case in 
which the kernel solves the Laplace's equation in the space of higher spatial di- 
mensionality, as is common in many other systems with competing interactions, 
see e.g. [7,16]). 
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Figure 1: A multi-droplet pattern: density plot of u in a local minimizer of 
£[u] with W(u) = j(l — u 2 ) 2 obtained numerically for u — —0.5, e — 0.025, 
and £1 = [0,11.5) x [0,10), with periodic boundary conditions. Dark regions 
correspond to u ~ —1, and light regions correspond to u ~ 1 (from [14]). 



The parameter e > in (1.1 1 determines both the scale of the short-range 
interaction and the magnitude of the interfacial energy between the regions 
with different values of u when e is sufficiently small. In fact, it is known 
that no patterns can form in the system if e is sufficiently large [13, 14, 38]. 
On the other hand, when e <C 1, the first term in the functional £ becomes 
a singular perturbation, giving rise to "domain structures" (see Fig. [TJ , which 
are of particular physical interest. These patterns consist of extended regions 
in which u is close to one of the minima of the potential W, separated by 
narrow domain walls. In this situation one can reduce the energy functional 
appearing in ( |1.1[ ) to an expression in terms of the interfaces alone. In [13,14], 
such a reduction was performed for £ using formal asymptotic techniques (see 
also [30,35,39,40]) and leads to the following reduced energy (for simplicity of 
notation, we choose the normalizations in such a way that the parameter e is, 
in fact, the domain wall energy, see Sec. [I] for details): 

E[u]= E - { \Vu\dx+\ ( [ (u{x)-u)G(x,y)(u(y)-u)dxdy. (1.4) 
1 Jo, z JnJn 

Here the function u takes on values ±1 throughout f2, and the kernel G is the 
screened Coulomb kernel, i.e., it solves the Neumann problem for 

-AG(x,y) + K 2 G(x,y) = 5(x-y), (1.5) 

with some n > 0. The constant k has the physical meaning of the inverse of the 
Debye screening length [13, 14]. Note that the sharp interface energy E with the 
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unscreened Coulomb kernel (i.e. with k = 0) was derived by Ren and Wei as the 
r-limit of the diffuse interface energy £ under assumptions of weak non-local 
coupling (i.e., with an extra factor of e in front of the Coulomb kernel) and 
u e (—1, 1) independent of e, as e — > [30] (see also [35]; note that this case is 
also equivalent to considering £ on the domain of size 0(e 1 ^ 3 )). At the same 
time, screening becomes important near the transition between the uniform and 
the patterned states which occurs near \u\ = 1, the case of interest in the present 
paper [13, 14]. Note that in the presence of screening the neutrality condition 



in (1.3 1 is relaxed. 

In this paper, we rigorously establish the relation between the sharp interface 
energy E and the diffuse interface energy £, and analyze the precise behavior 
of minimizers of the sharp interface energy E for e <C 1 in the vicinity of 
the transition from the trivial minimizer to patterned states occurring near 

= 1. We note that despite the apparent simplicity of the expression for E, 
the minimizers of E exhibit quite an intricate dependence on the parameters 
for e <C 1 and |u| ~ 1. Our analysis in this paper will be restricted to the 
case d = 2. While a number of our results can be extended to arbitrary space 
dimensions, our methods to obtain sharp estimates for the energy of minimizers 
rely critically on the properties of minimal curves in two dimensions and the 
logarithmic behavior of the Green's function of the two-dimensional Laplacian 
near the singularity. Therefore, they cannot be readily extended to other spatial 
dimensionalities, and, indeed, one would expect certain important differences 
between these cases and the case of two space dimensions. At the same time, 
we will show that in the case d = 2 it is possible to obtain rather detailed 
information about the structure of the transition near \u\ = 1 in terms of energy. 
Let us note that, since the case d = 1 is now well-understood [29-31,41], the 
remaining open case of physical interest is that of d = 3. 

Before turning to the analysis, let us briefly mention a perfect example of 
an experimental system in which the regimes studied by us could be easily 
realized, which is inspired by the beautiful Nobel Lecture of Prof. G. Ertl 
[42,43]. Consider molecules which undergo adsorption and desorption to and 
from a crystalline surface. On the surface, the atoms may hop around and 
reversibly stick to each other to form monolayer aggregates [44]. Then, within 
the framework of phase field models, this process may be described by the 
following evolution equation for the adsorbate density fraction <p [25] : 

(f> t = MA{W'{<j>) - gA(j)) + fc+(l -(f)- k-(f>, (1.6) 

where W is a double- well potential with two minima between </> = and <j> = 1, 
g is the short-range coupling constant, M is a kinetic coefficient, and k± are the 
adsorption and desorption rates, respectively. Note that this equation can be 
rewritten as 

4> t = MA{W'((f>) - gA</> + kG * (0 - $)}, (1.7) 
where k — (k + + fc_)/M, <j> = k + /(k + + £;_), and "*" denotes convolution in 



space, with Go given by (1.2|, provided the spatial average of the initial data is • 
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Upon suitable rescaling, this is precisely the H^ 1 gradient flow for the energy 
£, i.e., we have Ut = A(S£ /Su), where u is a rescaling of <fi. In particular, 
minimizcrs of £ are ground states of the considered system in equilibrium in 
the mean-field limit. We note that the adsorption and desorption rates k± can 
be quite small compared to the hopping rate, resulting in very small values of 
e ~ k 1 / 2 . Therefore, one can achieve a very good scale separation between the 
interfacial thickness (atomic scales) and the size of adsorbate clusters (micro- 
scale) in this experimental setup. 

Our paper is organized as follows. In Sec. [2] we present heuristic arguments 
and give the statements of main results, in Sec. [3j we perform a detailed analysis 
of the sharp interface energy E, in Sec. |4]we establish a connection between the 
sharp interface energy E and the diffuse interface energy £. Finally, in Sec. [5] 
we conclude the proofs of the theorems. 

Throughout the paper, the symbols W , H k , W k,p , C fc ' Q , BV denote the 
usual function spaces, | • | denotes the d-dimensional Lebesgue measure or the 
(d — l)-dimensional Hausdorff measure of a set, depending on the context, and 
C, c, etc., denote generic positive constants that can change from line to line. 
The symbols 0(1) and o(l) denote, as usual, uniformly bounded and uniformly 
small quantities, respectively, in the limit e — > 0, etc. Finally, we will say that a 
statement holds for e C 1, etc., if there exists e > such that that statement 
is true for all < e < e . For simplicity of notation, the subscript e is omitted 
for all quantities depending on s. 

2 Heuristics and main results 

Let us begin our investigation by setting d = 2 and making a simplifying assump- 
tion that the domain Q, is a torus: f2 = [0, l) 2 . Let us also specify the domains 
of definition for the functionals £ and E. Formally, the diffuse interface energy 
£[u) will be defined for all u € i? 1 (fi) subject to J n udx = u, whereas the sharp 
interface energy E[u) will be defined for all u G BV(il; {— 1, 1}). 

The assumption that 51 is a torus, which is common in the considered class 
of problems, eliminates the need to deal with the boundary effects and, even 
more importantly, restores the translational invariance inherent in the problem 
on the whole of Mr (note that the choice of the size of f2 is inconsequential, 
the obtained energy of the minimizers scales linearly with |f2|). As a result, the 
kernel of the non-local part of the energy becomes a function of x — y only. With 
a slight abuse of notation, in the following we will, therefore, replace G(x, y) 
with G(x — y) everywhere below. 

On heuristic grounds one would expect that the minimizers of E at e <C 1 
would be periodic with period R ~ e 1 / 3 , whenever |zt| < 1 and |u| is not too close 
to 1 [9,13,14,19]. A simple scaling analysis shows that in this case E ~ e 2 / 3 as 
e — > with u fixed. Our first result gives a justification for this energy scaling 
without any assumptions about the minimizers (for statements about existence 
and regularity of minimizers, see the following sections). 
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Theorem 2.1. Let W satisfy the assumptions (i)-(iv) at the beginning of Sec. 
[7| and let u € (—1, 1) be fixed. Then there exist e$ > and C > c > 0, such 
that 



cs 



2/3 < minE, min£ < Ce 2/3 



(2.1) 



/or aZ/ e < £o ■ 

Observe also that for E this result still holds when Q, — [0, l) d for any d, while 
for E it holds at least for d < 6 (see Sec. [4|. We note that for u — and \u\ < 1 
such a result was obtained by Choksi, using somewhat different techniques [9]. 
On the level of E (with k = 0), Alberti, Choksi and Otto recently proved, among 
many other interesting results, a stronger statement that in the limit e — > 0, the 



constants in the upper and lower bounds in (2.1 1 can be chosen to be arbitrarily 
close to each other [36]. We note that the case k = and u £ (—1, 1) fixed 
can be treated as the limit of energy E considered by us as k — > 0, when the 
constraint J^udx = u gets automatically enforced (see (5.2 1). 



Thus, when u 6 (—1, 1) is fixed, the energy E admits a non-trivial minimizer, 
whose energy scales as in (2.1 ) when e -C 1. What about the case |m| > 1? Here, 
in fact, it is easy to see that the only minimizcrs admitted by E are the trivial 
ones. Consider, for example, the case u < —1, the other case is equivalent 
by symmetry. In this case the problem admits the unique global minimizer 
u = — 1. To see this, let us introduce the characteristic function xn+ of the set 
tt + = {u = +1} for a given u G BV{£1; {-1, 1}). Then u = 2%n+ - 1, and by a 
straightforward computation 



E[u) > \ [ [ (2 Xn+ (x) 
z Jn Jn 



1 - u)G(x - y)(2xn+(y) ~ 1 ~ u)dxdy 



> 



(l + u) 2 _ 2(1 + 72) 
2k 2 k 2 



(2.2) 



Thus, when u < — 1, the second term in the last inequality in (2.2 1 is positive, 
hence, is minimized by |fi + | = 0. But in this case u — —1 attains equality in 
(2.2 I, so u = —1 is the minimizer. Thus, when |u| > 1, non-trivial minimizcrs 



of E do not exist, and, therefore, at \u\ — 1 a bifurcation occurs in the limit 
e -> 0. 

The main purpose of this paper is to investigate the transition between the 
trivial and the non-trivial minimizers of E and £ that occurs in the neighborhood 
of \u\ = 1 for £ <C 1. The energy E captures most of the difficulty associated 
with the considered problem. Therefore, we will spend most of our effort in this 
paper to the studies of E (see Sec. [3]). At the same time, as we show later 
(see Sec. [1]), the statements about the behavior of mini? also extend to that of 
min £ for £ < 1 (the correspondence of minimizers of the two energies will be a 
subject of future study) . 

When Q, — [0, l) 2 , the kernel G has an explicit representation 



G(x) 



nGZ 2 



(2.3) 
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where Kq is the modified Bessel function of the first kind. In particular, G > 
and we have the following asymptotic expansion from the power series repre- 
sentation of K [45] : 

G(x) = -^ln(R\x\) + 0(\x\), |x|«l, (2.4) 

where 



R=|/cexp 7 - J2 *o(«M) , (2.5) 

V n£Z2\{0} / 

and 7 « 0.5772 is the Euler's constant. We also have G(x) bounded whenever 



|x| > S, for any 5 > 0, and (2.4 1 can be used to estimate derivatives of G to 
0(\x\ In |x|) as well. 

Consider the case in which the value of u approaches u = — 1 from above, 
with £< 1 fixed. Clearly, for large enough deviations there exists a non-trivial 



minimizer. As can be seen from the arguments in the proof of Theorem 2.1 
the size of the set where u = 1 on the minimizer goes to zero as u —* —1. 
Heuristically, one would, therefore, expect that in this situation the minimizer 
will consist of a number of isolated droplets where u = +1 of small size in 
the background where u = — 1. Moreover, since on the scale of a droplet the 
interfacial energy will give a dominant contribution, these droplets are expected 
to be nearly circular. This motivates an introduction of the following reduced 
energy: 

JV 

E N ({n},{xi}) = ^(27rer J -2^(l + u) K - 2 r 4 2 - 7 rrf(lnKr,-|)) 

i=l 

iV-1 N 

+ 47f2 E E Gixi-x^rfrl (2.6) 

i— 1 

which describes the energy of interaction of N well separated disk-shaped droplets 
of radius r.j centered at Xj, to the leading order. More precisely, the first term 



(2.6 1 stands for the interfacial energy of all the droplets, the second term is the 
energy of interactions between the droplets and the background, the third term 
is the self-interaction energy of each droplet, and the last term is the interaction 
energy of each droplet pair (for the case of a single droplet in K 2 , see [14]). 



We can use the reduced energy in (2.6 1 to obtain the leading order scal- 



ing of various quantities for e <C 1 by balancing different terms. From the 
balance of interfacial energy and the self-interaction energy, one should have 
Ti = 0(e 1//3 | lne| -1 / 3 ). Balancing this with the second term leads, in turn, to 

6 = s- 2/3 \lns\- 1/3 {l + u) (2.7) 

being an O(l) quantity. Similarly, balancing the last term with the first three 
leads to N = 0(\ lne|), and the expected behavior of mini? at = 0(e 4 ^ 3 \ lne| 2 / 3 ). 
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One would also expect that, since the droplets repel each other, in a minimum 
energy configuration they would become uniformly distributed throughout f2. 

Our main result proves and further quantifies this heuristic picture on the 
level of the sharp interface energy E. 

Theorem 2.2. Let u — — 1 + e 2 / 3 | lnel 1 / 3 ^, with some 6 > fixed. Then for 
any a > sufficiently small there exists £q > such that for all e < e$: 

(i) If 5 < ^\/9n 2 , then u = —1 is the unique global minimizer of E, with 
e- 4 / 3 |lne|- 2 / 3 min£;= \k~ 2 P . 

(ii) If 6 > there exists a non-trivial minimizer of E. The minimizer 

is 

N 

u(x)=-l + 2j2xn+(*)> ( 2 - 8 ) 

8=1 

where Xfi+ ire characteristic functions of N disjoint simply connected sets 

flf C fl with boundary of class C 3 , and N = 0(|lne|). The boundary of 
each set flf is 0(e 2 / 3 ~ a )-close (in Hausdorff sense) to a circle of radius 
Ti centered at Xi. Furthermore, 

mhiE= ±e 4/3 \ln e\ 2/3 k- 2 6 2 + EN^n},^}) + 0(e 5/3 - a ), (2.9) 

with E N = 0(e 4 / 3 |lne| 2 / 3 ) ; n = O^^lne]' 1 / 3 ) , and 

\x i -x j \>e a , Vj^i. (2.10) 

(Hi) If 5 > | v9 k 2 , in the limit e —> we have 

£- 1/3 |lne| 1/3 r 4 -> (2.11) 

uniformly, 

N 



weakly in the sense of measures, and 




e- 4 / 3 |ln £ r 2 / 3 min£ - — I 6 - — k 2 ) . (2.13) 



Note that a more detailed result on the s tructu re of the transition occurring 
near § = i \/9 k 2 is presented in Proposition 



3.20 



Let us make a few remarks related to the statements of Theorem 12.21 For 
small, but finite values of e this theorem establishes an equivalence between the 
sharp interface energy E and the energy of N interacting droplets En, in the 
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N =2 



N =3 



N =4 



N = 5 



Figure 2: "Coulombic dice": Minimizers of E N with r t = v/3 e 1 / 3 1 In e| 1 / 3 for 
K = 2 and N = 2, 3, 4, 5, obtained using random search algorithm. 



sense that the minimizers of E are close to "almost" minimizers of En, i.e., we 
have Ejq < mini^ + 0(e 5 / 3_cr ). Nevertheless, to prove closeness of minimizers 
of E to those of En we also need some coercivity of the energy En- This problem 
has to do with the properties of the minimizers of the pairwise interaction of the 
droplets, i.e. the choice of Xi which minimize En with fixed r». This becomes 
a difficult problem in the case of interest, since we generally expect N 3> 1 (for 
a numerical solution at a few values of N and n — 2 see Fig. [2]). It would be 
natural to conjecture that at small enough e the minimizing droplets will arrange 
themselves into a periodic lattice close to a hexagonal (close-packed) lattice. 
Proving this kind of result, however, is a major challenge (see [46] for a recent 
proof for a certain class of pair interactions), which is one of the open questions 
also in many other problems, such as the problem of characterizing Abrikosov 
vortex lattices, for example [47]. Let us mention here a recent result by Chen 
and Oshita, who proved that in the case k = the hexagonal arrangement of 
disks is energetically the best among simple periodic lattices [48]. Yet, it is not 
known if the same result also holds for more general arrangements of droplets. 
Here we prove a weaker result that the number density of droplets becomes 
asymptotically uniform as e — > 0, leading also to uniform distribution of energy 
(compare with [36]). Moreover, we identify the precise asymptotic behavior of 
the minimal energy and the size of the minimizing droplets as e — > 0. 

Lastly, we establish the asymptotic behavior of the minimal value of the 
diffuse interface energy £ (which, of course, agrees with the result for the sharp 
interface energy). 

Theorem 2.3. Let W satisfy assumptions (i)-(iv) at the beginning of Sec. [7 ) 
let u = — 1 + e 2 / 3 | Ins] 1 / 3 !, with some 5 > fixed, and let k be given by (4-10). 
Then 

(i) IfS< \?/§k 2 , then £- 4 / 3 |ln£|- 2 / 3 min£ -> \kT 2 P , 

(ii) If~5>\^K 2 , i/jen£-V3|ln £ |- 2 / 3 min£ - 4^ 
as £ — > 0. 



The proofs of Theorems |2.1 2.3 are based on a number of propositions es- 



tablished in Sees. [3] and [4] and are completed in Sec. [5] 
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3 Analysis of the sharp interface problem 

Our plan for the analysis of the sharp interface problem consists of a number of 
steps which we list below: 

1. Introduce a suitably rescaled energy E and domain Cl. 

2. Establish existence and regularity of the minimizers of E (subsets of f2 
where u = 1). 

3. Establish some a priori estimates for the geometry of the minimizers of E 
and uniform bounds on the induced long-range potential. 

4. Establish that different connected components of minimizers of E are sep- 
arated by large distances in f2. 

5. Establish that each connected component of a minimizer of E is close to 
a disk (hence the term "droplet" ) . 

6. Establish equivalence between min E and min E^ (the suitably rescaled 
version of En). 

7. Improve the estimate for the separation distance between different droplets. 

8. Prove uniform convergence of the rescaled droplet radii to a universal 
constant. 

9. Prove convergence of min £ to a limit and convergence of the normalized 
droplet density in the original, unsealed domain 51 to a limit, as e — > 0. 

This plan is carried out in the rest of this section via a series of lemmas and 
propositions. 



3.1 Scaling 

We begin by introducing a suitable rescaling, in which the main quantities of 
interest become O(l) quantities in the limit e — > 0. Motivated by the discussion 
of Sec. [2] we define the rescaled energy E (with the energy of the uniform state 
u = —1 subtracted) and a new coordinate x € fl — [0, e 1 / 3 1 lne| -1 / 3 ) 2 , where Cl 
is a two-dimensional torus with period e 1 / 3 ! lne| -1 / 3 : 

E = £ 4 / 3 | lne| 2 / 3 (^P + E), x= x. (3.1) 

The energy E can be conveniently expressed in term of the set Cl + C Q in which 
u= 1: 

E = llnel^flan+l -2Sn- 2 \n+\) 

+ 2\lne\- 2 f f G(e 1/3 |lner 1/3 (x -y))dxdy. (3.2) 
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We also need an expression for the rescaled energy En of a system of in- 



teracting droplets. With the help of (3.1|, we can write the rescaling of (2.6 1 
as 

9 N 

' ' 2 — 1 

2 N ~ 1 N 

+ t^EE^ih- 1 ^-^ (3-3) 

where fj and Xi are the radii and the centers of the droplets, respectively. 
3.2 Properties of minimizers 

Let us begin with the statement of a result on the existence and regularity of 
minimizers of E (or, equivalently, of E), which is obtained by straightforwardly 
adapting the results of [49] for sets of prescribed mean curvature. 

Proposition 3.1. There exists a set f2 + of finite perimeter which minimizes 
E in (3.2). The boundary dCl + of this set is a curve of class C 1,a for some 



a e (0,1). 

In view of this, in the following we will always assume that minimizers f2 + of E 
are closed sets. We also note that 

v(x) = \lne\- 1 f G{e 1/3 \h-ie\- 1/3 (x-y))dy (3.4) 

is in W 2 ' p (Cl), with any p > 1, and, hence, in C 1,Q (f2) for any a E (0, 1). Indeed, 
v solves the equation 

- Av + e 2/3 \lne\- 2/3 K 2 v= |lne|- 1 xn+, (3-5) 

where Xn+ ^ s the characteristic function of f2 + , in f2, and so the result follows 
by standard elliptic regularity theory [50]. As a consequence, we have a higher 
regularity for the boundary of the minimizer f2 + of E ( [51], see also [35]): 

Corollary 3.2. The boundary d£l + of a minimizer Cl + of E is of class C 3,a . 

Note that this regularity result also holds more generally for local minimizers of 
E in dimensions d < 7 [49], hence, in particular, the expressions for the first and 
second variation of E in d < 3 obtained in [14] are justified (see also [52] for the 
case of arbitrary dimensions). If p € C 1 (dCt), a > 0, and Cl a is the set obtained 
by displacing dCl + by ap in the outward normal direction, then a i— > E(Q,j) is 
twice continuously differentiable at a = 0, and we have [14] (for the reader's 
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convenience, the computation is reproduced in Appendix O 



me 



dE(ti+) 



In el 



da 



(K(x) - 25k- 2 + Av(x))p(x) dH 1 ^), (3.6) 



da 2 



3Q+ 



(\Vp(x)\ 2 + Av(x) ■ Vv(x) p 2 (x)) dn\x) 



^llner 1 



an+ Jan+ 



G(e 1 / 3 |lne|- 1 / 3 (5-y))p(x)p(y) ( iH 1 (i) ( iH 1 (y). (3.7) 



where if (2) is the curvature at point x S <9f2 + , with the sign convention that 
if > if f2 + is convex, and is the outward unit normal to dti + at that 
point. The associated Euler-Lagrange equation for dti + reads 



K(x) = 2Sk~ 2 -Av(x), 



(3.8) 



which also allows to simplify the expression in (3.7 1 evaluated on a minimizcr 
to 

d 2 E(Cl+) 

a=0 



In el 



da 2 



{\Vp(x)\ 2 + Av(x) ■ Vv(x) p 2 {x) - K 2 (x)p 2 (x)) dn\x) 

dn+ 

+4|lne|- 1 / / G(e 1 / 3 |lner 1 / 3 (x-y))p(x)p(y)dH 1 (x)dH 1 (y). (3.9) 
Jan+ Jdn+ 

We will use these equations later on to establish some properties of the 
minimizers for e <C 1. Meanwhile, let us begin our analysis with some basic 
estimates. 

Lemma 3.3. Let fl + be a minimizer of E. Then there exists C > such that 

< C|lne|, (3.10) 
\dQ+\ < C\lne\. (3.11) 

for £«1. 



Proof. First of all, by representation (2.3 1 we have G(x — y) > c > for all 
x,y € CI. Therefore, in view of the fact that min E < (since E = if Ct + = 0), 
from (3.2 1 we have 



0>\1ne\E 2|lne 



-l 



iq+ Jn+ 
2c|lne|- 1 \Cl+\ 2 -2Sk- 2 \Ci + 



G(e 1/3 | lner 1/3 (x - y))dx - 2fe- 2 |fi+| 

(3.12) 



which gives (3.101. On the other hand, we also have 

< 28n- 2 \n+\. 



Therefore, from (3.10) we immediately obtain (3.111. 



(3.13) 

□ 
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As a corollary, it follows from ( |3.11 I that the diameter of each connected 
subset flf of f2 + is bounded by 0( \ In e| ) 



diam(fit) < C|lne|, 



(3.14) 



for some C > independent of e <C 1. 

Our next step is to show that the area of each connected component of 
Cl + =/= is uniformly bounded above and below independently of e. 

Lemma 3.4. Let Cl+ = u£L]fi+ be a non-trivial minimizer of E, where £lf are 
the disjoint connected components ofCl + . Then, there exist C > c > such that 



< |Q+|, |<9Q+| < C, diam(Q+) < C, 



(3.15) 



for e <C 1 . 



Proof. First, note that sin ce by Corollary 3.2 the set dil + is of class C 3,a we 
have N < oo. To see that (3.151 holds, we first write E as 

■2SK- 2 \nf\ 



\\ne\E = 




21 In el 



_1 E 



G(e 1/3 |lne|- 1 / 3 (x - y)) dxdy 

G(e 1/3 | lne|- 1 / 3 (5 - y)) dxdy ) . (3.16) 



In view of (3.14) and ( |2.4| , the integral in the second line in (3.161 is bounded 
from below by g^(l — S)\ lne| \Qf \ 2 for any S > 0, provided e is small enough. 
Therefore, removing the set flf from £l + will result in the change of energy AE 
estimated as 



| In el AE 



< 



\dn+\ -2~8k~ 2 pf\ 4 
2^m+| 1 / 2 -2fc- 2 i" 



+ 3^(1 -mt? 



(3.17) 



where in the first line we took into account that G > and in the second line 
used the isoperimetric inequality. Then, by direct inspection (see also Fig. |3| 
we have AE < 0, contradicting minimality of E on f2 + , unless c < |f2 + | < C for 
some C > c > 0, independently of e <C 1. Finally, the lower bound for |<9^ + | 
follows from the isoperimetric inequality, and the upper bound is obtained by 
applying the previous argument to the first line in (3.17). □ 

Following the same arguments, we also immediately arrive at the following 
non-existence result: 

Proposition 3.5. Let 5 < ^\/9n 2 be fixed. Then the unique minimizer of E 
is 0+ = for e < 1. 
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Proof. Let us introduce the function V v : [0, oo) — > R, defined as 

T4(f) = 2vr (f + (2w - (5K" 2 )r 2 + |r 4 ) , (3.18) 

whose graph at v — and several values of 5 is shown in Fig. [3] If f2^~ is a 
connected component of f2 + and fj = (^Ifj^l) 1 / 2 , then by the same arguments 
as in Lemma 3.4 the energy gained by removing Clf from fl + is bounded below 
by I lne| _1 (Vo(r.() + o(l)), as long as e < 1. Then, by direct inspection Vo(f) 
is always positive under the assumptions of the proposition, making 17 + = 
energetically preferred. □ 



Note that the asymptotic value of the threshold of 5 in Proposition 3.5 below 
which no non-trivial minimizers are present was computed in [14]. Another 
simple corollary to Proposition |3.4| is the following 

Lemma 3.6. Let Cl + be a non-trivial minimizer of E and let N be the number 
of disjoint connected components o/fi + . Then there exists C > such that 

7V<C|lne|, (3.19) 

for £«1. 

Let us now establish a uniform bound on the potential v. Note that a version 
of this result is also an important component in the proofs of [36] . 

Lemma 3.7. Let Cl + be a non-trivial minimizer of E. Then for any a € (0, 1) 
we have 

0<v<C, |Mlci.«(fi)<C (3.20) 



where v is given by (3.4-), for some C > independent of e <C 1. 
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Proof. We start by noting that v > in view of positivity of G. Let us now 



estimate the gradient of v. Using (2.4 1 and Lemmas 3.4 and 3.6 we get 



\Vv{x)\ < liner 1 / |VG(e 1/3 |lne|- 1 / 3 |a;-y|)|dy 
Jn+ 

1 / \VG(eV 3 \]ne\-V 3 \x-y\)\dy 

JB f (x) 



< |lne 
liner 1 



|VG(e 1 / 3 |lne|- 1 / 3 |x-y|)|dy 

n+\B f (x) 

< CQlne^r + f- 1 ) < 2C| lne^ 1 / 2 , 



(3.21) 



for some C > 0, where B f (x) is a disk of radius f centered at x, and the last 
inequality is obtained by choosing f = | In e| x / 2 . Therefore, by the results of 
Lemma 13.41 we see that 



osc u(x) 



(3.22) 



for each connected component (if of (l + . To see that this implies the conclusion 
of the lemma, suppose that, to the contrary, we have maxu = M 3> 1. Since 
by (3.5 I the function v is subharmonic in fA£2+ it achieves its maximum in the 
closure of some (if. Therefore, in view of (3.22) we have v > \M in (if. Then, 
following the same arguments as in the proof of Lemma |3.4[ for large enough 
M we can lower the energy by removing (if from (l + . 

Finally by [50, Theorem 9.11] we have |?.> | w 2 -p{b 1 (s)) ^ C; where Bi(x) is a 
disk of radius 1 centered at x € Q, for some C > and any p > 2, independently 
of x and Hence, the uniform Holder estimate on the gradient follows by 

Sobolev imbedding. □ 



We can also immediately conclude from (3.8) and (3.221 that the curvature 



of dfl + is uniformly bounded both from above and below by positive constants, 
implying that each flf is convex. Note that this result justifies the terminology 
"droplet" for each Clf which we will be using from now on. 

Lemma 3.8. Let dfl + be the boundary of a minimizer fl + of E. Then we have 

c < K(x) < C, (3.23) 

for all x € d£l + , with some C > c > independent of e <C 1. In particular, when 
each connected component (if of(l + is convex and simply connected. 



Proof. The upper bound is an immediate consequence of (3.8 1 and positivity of 



v. To obtain the lower bound, let us note that by the results of Lemma [3T4[ for 
every connected component Clf there exists a disk Bf.(xi), with ?i = O(l), such 
that flf C B fi (xi). Therefore, translating _B fi (xi) until its boundary touches 
dflf , we obtain a point x\ € dfl f , s uch that -ftT(x^) > r" 1 > 2c, for some c > 
independent of e 1. Now, by (3.8 1 we have u(x£) < \{5n~ 2 — c). At the same 



time, by (3.22) this implies that v(x) < \(2Sk 2 — c) for all x € dClf, which, 
again, by (3.8 1 gives the statement. □ 
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We now show that different connected components of f2 + cannot come too 
close to each other when e <C 1. 

Lemma 3.9. Let Cl+ = U^Clf be a non-trivial minimizer of E, where £1^ 
are the disjoint connected components ofCl + , and let N > 2. Then, there exists 
C > such that 

distant) > C Vijt j, (3.24) 

for £«1. 

Proof. Let x t £ tif and Xj € be such that r = — Xj \ = dist (f2^~ , f2+) > 0. 
Consider a disk B centered at \{x\ + X2) with radius R — 2r and a rectangle 
Q inscribed into B which is shown by the thick solid lines in Fig. [4} In view of 



the uniform bound on the curvature of <9£! + obtained in Lemma 3.8 the curve 
segments d£lfr\Q and d£l~jT\Q passing through and Xj, respectively, intersect 
dQ transversally as in Fig. [4] when r 1. Furthermore, we have dist(<9f2^~ n 
dQ+,dn+ndQ+) < 2r and dist(dttf ndQ~ , dCtf ndQ~) < 2r, where dQ+ and 
dQ~ are the right and the left side of the boundary of the rectangle relative to 
the line through x\ and £2, respectively, for sufficiently small r independent of 
£ < 1 (see Fig. At the same time, we have \dflf (1Q\ + \dfl+ (1Q\> 4r%/3. 
Therefore, reconnecting the points d£lfr\dQ + with dtij C\dQ + , and dClf(ldQ~ 
with dCt^ (~l dQ~ by straight lines and including the region between them into 
f2, we will decrease |<9f2 + | by at least 4(v / 3 — l)r. Thus, the change AE in the 
total energy is estimated to be 



|lne|A£ < -4(V3- l)r + 4 / v(x)dx + 

Jo 



2|lne|- 1 / / G(e 1/3 \lne\- 1/3 (x-y))dxdy. (3.25) 
Jq Jq 



Finally, in view of Lemma 3.7 and (2.4 1, the the right-hand side of (3.251 is 



bounded above by — C\r + Cir 1 , with > independent of e <C 1. Hence, 
the energy of such a rearrangement will be lower if r is sufficiently small, for all 
e < 1. □ 

As our next step, we establish that different droplets must, in fact, be suffi- 
ciently far from each other. We note that this result is a manifestation of the 
"repumping" instability, which does not allow two droplets to approach each 
other sufficiently closely. Dynamically, this instability results in the growth of 
one droplet at the expense of its neighbor shrinking. This instability mechanism 
for reaction-diffusion systems was first pointed out in [53] (see also [4]) and fur- 
ther studied in the context of two-dimensional periodic structures in [13,14,54]. 

Lemma 3.10. Let n+ = uf =i nf be a non-trivial minimizer of E, where 
are the disjoint connected components of £l + , and let N > 2. Then there exists 
a > such that 

dist(fit, n+) > £ - a Vi + j, (3.26) 



16 




Figure 4: Schematics of the rearrangement argument of Lemma 3.9 In (a), the 



set f2 + is shown in gray, solid arcs show the bounds on the location of <9f2 + , 
the thick solid lines show the rectangle Q. In (b), the gray region shows the 
rearranged Cl + . 



for e <C 1 . 



Proof. Consider the second variation of E with respect to the perturbation, in 
which the boundary of each connected component tif is expanded uniformly by 
a distance aci in the normal direction, i.e., we have p(x) — Ci for all x £ dft^ . 
By (3.9 1, we have 



d 2 mt) 



da 2 



= \ hxe\ 1 ^QijCjCj, 



(3.27) 



where the coefficients Qij of the quadratic form Q can be estimated as 



K 2 (x)dH l (x) + — \dn 
+ Sir 



+ 12 



o(l), 



(3.28) 
(3.29) 



where we took into account that by (3.5 1 and Gauss's theorem Jqq+ v{x) ■ 

Vv(x) dH 1 (x) — — | lnel -1 !^^ | + 0(e 2/3 | lne| -2 / 3 ) and used the expansion in 
(2.4 1 together with Lemmas 3^8 |3.7 and 3.4 for e«l. Furthermore, since by 
convexity of flf (see Lemma 3.8 1 the boundary of each flf is a closed curve, by 



Cauchy-Schwarz inequality we have 



4tt 2 = 



KffldH 1 ^)) <\dn+\ K 2 (3:)dM 1 (x). 



(3.30) 



Therefore, the diagonal elements of Q can be further estimated as 



on 



4tt 2 



0(1). 



(3.31) 



On the other hand, define a = |lne| 1 ln(dist(fi^, Qf)), and suppose, to the 
contrary of the statement of the proposition, that a is sufficiently small for some 
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pair of indices for a sequence of e — ► 0. Then, with the help of Lemma 3.9 and 

2 



(2.4 1 we can estimate 



Q iJ = — (l-3a)|flfi+||an+|+o(l). 



3tt 

Now, for the index pair above let us choose c, = |<3ffi 



(3.32) 

,- = -WI, 

and let us set — for all other indices k. A simple calculation of the sum in 
(3.271 then shows that for this choice of c's we have 



hi; 



da 2 



< 



A\dn+\ 2 \dn+\ 2 



o(l), (3.33) 



where we took into account Lemma 3.4 This expression is negative for e -C 1, 
if 



a<2-K s mm{\d9+\- i ,\d9+\- i }, 



(3.34) 

which, in view of Lemma |3.4| contradicts minimality of E for small enough 
e. □ 



Let us also point out that the proof of Proposition 3.10 gives a universal lower 
bound for the perimeter of the connected portions of the minimizers. Indeed, 
the quadratic form Q has a negative eigenvalue, if ^\d£lf | 2 — A-K 2 \d£lf \~ 1 < 
and e -C 1, which implies the following result: 

Proposition 3.11. Let Cl+ = U^ =1 Cl+ be a non-trivial minimizer of E, where 
flf are the disjoint connected components of fl + . Then, for every 8 > 

\dn+\>TT^6-6, (3.35) 

for £«1. 

Note that this criterion in the radially-symmetric case was obtained in [13,14,55] 
and is also applicable to all local minimizers (for global minimizers, a better 
bound will be obtained below) . We also derive another quantitative estimate on 
v and the geometry of Clf that remains valid for local minimizers of low energy. 

Proposition 3.12. Let 0+ = uf =i nf be a non-trivial minimizer of E, where 
Clf are the disjoint connected components of Ct + . Then 



< v < 



2k 2 : 



in? 



< 



3tt8 



(3.36) 



for £«1. 

Proof. Let x £ dClf. Then, using Lemma 
for some c > independent of e <C 1 we obtain 

< c < K (x) 



(3.81, (2.4), and positivity of G, 



< 25k 



41 lnel 



2Sk~ 2 - Av(x) 
G^llner^i-y))^ 
2 



< 26k~ 



3tt 1 



(3.37) 
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which, together with (3.22) and the fact that v reaches its maximum in the 
closure of f2 + , yields the statement. □ 

We now prove that for e<l each droplet in a minimizer is, in fact, close to 
a disk. The basic idea of the proof is that because of the logarithmic behavior 
of G at small distances the potential v inside each droplet is approximately 
constant. Therefore, the shape of the droplet approximately minimizes the 
usual isopcrimctric problem, and the size of the droplet is determined by the 
balance of surface tension and the pressure due to non-local forces inside the 
droplet [13,14]. 

If the droplet £lf were exactly a disk Bf^Xi) of radius fi centered at £j, 
then the potential v would be given by 



v*(x) = v*(x) + Vi(x), 



where 



nSZ 2 



v B (\x - Xi 



with the function v B (p,r,K) being the solution of — Av B + k 2 v b = 
K 2 , given explicitly in terms of the modified Bessel functions: 



v B (p,r, k) 



K — k rKi(Kr)Ia(Kp), p < r, 



k 1 rI 1 (Kr)K (np), 



and 



v i (x) = \lne\- 1 J2 [ G^lne^ix-y))^. 



■ / ■ J (I ■ 



Note that in view of Lemmas 3.10 and 3.4 and (2.4), we have 



IVwil < Ce a , 



(3.38) 

(3.39) 

Xs r (o) in 

(3.40) 
(3.41) 
(3.42) 



for some C > and a > 0, in any disk of 0(1) radius containing flf for 
Therefore, if x £ dBf^Xi), by Taylor-expanding the Bessel functions [45] we 
have for any a < \ 

v*(x)=v l + 0(e a ), |V«*(2)| ^Odlner 1 ), (3.43) 

where the constant Vi is given by 

v 4 = - illnel-iffln^llner 1 / 3 ^) 

+n\lns\- 1 Y, fp(e 1/3 nns\- 1 ^ K (x l -x J )) : f 3 - = (^t\) 1/2 . (3.44) 

In the following, we will show that v(x) on dflf also coincides with Vi to 0(e a ), 
giving the balance of forces at the interface. We are now ready to state our 
result: 
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Proposition 3.13. Let n + = U^fit be a non-trivial minimizer of E, where 
VL~i are the disjoint connected components of f2 + . Then there exists a constant 
a > such that for all e <C 1 .' 

(i) For each (if there exists a point Xi € (if, such that 

B fl - e »(xO C n+ c B f . +ea {xi), (3.45) 

where fj = (^If^l) 1 / 2 and B f (x) is a disk of radius f centered at x; 

(ii) The values of fi satisfy 



25h 



(3.46) 



where Vi are given by (3. 44) 



Proof. Let us pick a point x\ S fi^ , then 0+ C -B| a ^+| (a;-)- Let us then replace 



Qf with a disk of the same area centered at x\. By Lemmas 3.10 and 3.4 the 
resulting set B fi {x') still satisfies the bound in (3.261, and the change of energy 



AE under this rearrangement can be estimated as 



|lne|A£; = 2^1^: 



|l/2 



|dQ+| + 0(|lm 



(3.47) 



where we used (2.4 1, (3.421 and Lemma 3.4 Thus, the energy will decrease 



under this rearrangement, contradicting minimality of E, unless for some C > 
the isoperimetric deficit of (if 



i>(n+) = 



\dSl1 



- 1 < C|lne| 



(3.48) 



for e <C 1. Choosing X{ £ ^ign+l ( x 'i) t° minimize \£lfABr i (xi)\, where £lf AB fi (xi) 
denotes the symmetric difference of sets (if and S fi (5j), by the results of 
[56] we have \(lf AB f .(xi)\ < C'\ In el -1 / 2 , and C > is a constant indc 



pendent of e -C 1. Then, by Lemma 3.8 the set dflf is uniformly close 
to dBf^Xi), and 5$ € (if (if not, then by convexity of (if we would have 
\n+ABr z ( Xl )\ > l\B fi ( Xi )\), giving (i) to o(l). 

To obtain the 0(e a ) bound in (i), let p : dB fi (xi) — > K be the signed 
distance from a given point on dB f .(x.i) to <9f2^~ along the outward normal to 
dBf.(xi). Note that by convexity of (if the function p defines a one-to-one 
map between dflf and dB fi {xi). Furthermore, if ||/9||i«>(aB ? .(s j )) = 6, we have 
11^ ' p\\L°°{dBr (xi)) < C5 1 / 2 for some C > and e -C 1 in view of Lemma 
► 0, as £ — > 0. Also, by Corollary 3.2 we have p € C 3 (<9i?r; (a^))- 



and 5 



3.8 



Treating fl + as a perturbation of the set fl* = B fi (xi) U (f2 + \f2^) and 
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expanding as in Lemma C.l we can write 



]ne\(E(Q+) - E(ST)) 



dB n {Si) 



(f" 1 - 25k~ 2 + iv*(x)) p(x) dTi}{x) 

dB f . (Si) 

(|Vp(i)| 2 + Av{x) ■ V/(i) p 2 (x)) dH 1 ^) 



1 

'2fi 



i JdB f .(xi) 



(4v*(x) - 2Sk~ 2 ) p 2 (x) dH\x) 



2 



OBf.(xi) JdB f .{xi) 



G(e 1 / 3 | lnsl" 1 / 3 ^ - y) p(x)p(y) dH^dH^y) 



0(^ +Q ),(3.49) 



for any a S (0, 1). Moreover, in view of Lemmas 3.7 and 3.8 the error term in 
(3.491 is uniform in e <C 1. 

On the other hand, since |f2^~| = \B fi {xi)\, we have 







dBf.(xi) JO 



p{x) 



(1 + f t "V) drdH l {x) 



dB?.(xi) 



pi^dH 1 ^) 



2r, 



i JdBf.(xi) 



p 2 (x)dH 1 (x). (3.50) 



Therefore, using the estimate in (3.431 we can rewrite (3.491 as 



| lne\(E(Q + ) - E(Q*)) 



1 



dBf.(xi) 



(\Vp\*-rrVW)dH\z) 



\lne\ 



B Fi (Xj) JdBf^Si) 

0(| ^e\~ l \\p\\h 



G(e 1 / 3 | lnel" 1 / 3 ^ - y) p(x)p(y) dH^^dU 1 ^) 



0(s a \\p\\ Ll _ „ J+0(S a \\ P \\ 2 Hi ,_ J,(3.51) 



where we took into account that 8 < C\\p\\ H i for some C > 0. Further 
estimating the double integral in ( 3.51| , using 



C<51n|lne| 



we have 



In el 



|lne|(£(ft+) -£(«*)) 



(3.52) 



as f . (xi) 



{\Vp\ 2 -r7 2 p 2 (x))dH 1 (x) + 



3tt 



SBp. (Si) 



p&dH 1 ^) 



-0{e a \ 



2 



HI r .)■ (3-53) 
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Now, write p as p = p + pi + P2, where p = JoB^(xi) P( x ) c ^ 1 ( a; ): 
Pi(x) = rfEf 1 ] ' b, for some vector 6 e I 2 , and P2 orthogonal to po and p\ in 



L 2 (dB f% {xi)). By (3.501 we have \p \ = 0{\\p\\ 2 L2{dBf (s . z)) ), which is, therefore, 
negligibly small compared to |6| and ||p2||L 2 (9B F .(2 i )) in all the arguments below. 
Then, using Poincare's inequality, we find that 



E(Q + ) > E(Q*) 



\Vp 2 \ 2 {x)dH 1 {x)-Ce a \\p\\ Ll 



IdBf. (Si) 

for some C > and < c« 1, whenever e < 1. This implies that 
INrk <4C^||p|U +4c|&| 2 , 



(3.54) 



(3.55) 



for e <C 1, otherwise replacing £lf with B f .(xi) lowers the energy. On the other 
hand, we also have |fe| = OGI/^Hh 1 )■ If not, then dClf will be o(\b\) close to 

dB fi (xi + b) for e<1. This, however, contradicts the choice of Xj to minimize 
\9,fAB f .(Xi)\. Therefore, we have 



Hi 



< Ce a \ 



c\\p\\hI 



(3.56) 



for some C > and < c < 1, implying ||p||#i = 0(e a ) and, hence, 

5 = 0(e a ). This gives part (i) of the statement of the proposition. 

Finally, to prove part (ii) of the statement, let f2+ be obtained from Cl + by 
expanding flf by an amount a > 0, i.e., let us change p(x) — > p(x) + a for every 
x S dB fi (xi). By (3.49), the change of energy can be estimated as 

\lne\(E(n+) - E(n+)) = l^af^f' 1 ~ 2<5«T 2 + Av t ) + 0(aS) + 0{a 2 ), (3.57) 



where we took into account (3.43). Then, since £l + is a minimizer, the right- 
hand side of (3.57) should vanish to O(a). Therefore, by previous result we 



obtain the statement. 



□ 



Also, from the proof of Proposition |3.13| we obtain the following universal 
lower bound on 

Proposition 3.14. Let n + = U^fit be a non-trivial minimizer of E, where 
Clf are the disjoint connected components of fl + . Then, for every S > 



(3.58) 



for £«1. 



Proof. By Proposition 3.13 each value of = (^jf^l) 1 / 2 satisfies (3.461 and, 



hence , is cl ose to a critical point of V v jxt) defined in (3.181, which can be seen 
from (3.44) by Taylor expansion. Furthermore, we must have 2vi(xi) — 5k~ 2 < 
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— \\f§ + o(l), so that V Vi (xi)( r i) < f° r £ < 1| otherwise, arguing as in 
Lemma 3.4 we can reduce the energy by removing 0,f from f2 + . Therefore, 
fj should be close to the positive minimum of V Vi ( Si y By inspection, in this 
situation f, > \/% — 5, for any 6 > (see Fig. [3]), provided e is sufficiently small, 
hence, the claim. □ 



The results of Proposition 3.13 just obtained immediately allow to establish 
an asymptotic equivalence of the energy E and the reduced energy En on the 
minimizers for 



Proposition 3.15. Let VI 

Proposition \3l3\ Th 

minis = 0(1), 



+ 



be a non-trivial minimizer of E, where 
are the disjoint connected components of Cl + , and let fj and Xi be as in 
en 



min E = min E 



N 



0(e a ), 



(3.59) 



for some a > independent of e <^ 1. 



Proof. Th e first equation in (|3.59| is a direct consequence of the definition of 

, ■■ ■ ■ - 4/31 



E in (3.2 1, according to which < 



mini; < £- 4 / 3 |lne| 2 / 3 £;[-l] 



The upper bound in the second equation follows by choosing a trial 



function for E in the form of disks of radius r. 
En and taking into consi derati on Lemmas 3.4 



other hand, by Proposition 3.13 i), we have Ci^ 
\dQf\ > 2 7 r(r l 



centered at x 
and 



which minimize 
3.101 and (Ell. On the 



D B fi _ ea (xi) for s <C 1, hence, 



e Q ) and |f^| > n(fi — e a ) 2 . This controls from below all the 
terms of E, except the one involving 5, by the corresponding terms of En- The 
latter, however, is controlled by the second inclusion in Proposition 3.13 i) 



□ 



To summarize, for < e <C 1 the non-trivial minimizers of E have the form 
of well-separated nearly circular droplets. In fact, from Proposition |3.15| one 
should expect that the droplet-droplet inter actio n part of the energy, which 
is given by the last term in the expression (3.3) for En, should be close to 
the minimum for fixed droplet sizes. Proving this, however, generally requires 
information about coercivity of the interaction energy, which becomes difficult 
to establish when N ^> 1, the asymptotic case of interest. Nevertheless, with 



the help of Lemma 3.7 we can prove that in the original scaling the droplets 
stay away from each other a distance O( e^) i n fi, with an arbitrary f3 > for 
e <C 1, i.e., that the statement of Lemma 3.10 actually holds for any a £ (0, |), 
provided that e is small enough. 



Proposition 3.16. Let tt+ = U^fi^ be a non-trivial minimizer of E, where 



12+ 



3.13 



are the disjoint connected components offl + , 
Then, for any a G (0, i) we have \xi — Xj\ 



and let Xi be as in Proposition 
> £~ a , for all i j, as long 



as e <C 1. 



Proof. First of all, note that by Lemma 3.10 the statement of the Proposition 
holds for some a > 0. To prove that a could be chosen arbitrarily close to |, 
suppose that, to the contrary, there exists a sequence of e — > and a pair of 
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indices depending on e, such that Xj\ < e~ a with some < a < |. Let 
us denote by I\ the set of indices of those droplets whose centers are contained 
in a disk B\ centered at xo — \(xi + Xj) with radius e~ a . By assumption we 
have |Ji| > 2, where | • | denotes the counting measure. Also, we have < M 
for some M e N independent of e <C 1. Indeed, by Lemmas |3.4| and |3.7[ and by 



(2.4 1 we have for some c > 

C^^SoJ^llnel- 1 V / Gis^lhxel-^ixi-y^dy 



fce/i ' 



> (|-a + o(l))c|7i|, (3.60) 



for £ < 1. 

Now, fix c > sufficiently small independently of e, and consider a sequence 
of nested disks Bk of radii £- a ( 1 + ka ) centered at Xq. By repeating the argument 
above, we also have \Im\ < M, as long as e <C 1, where is the counting 
measure of the set Ik of indices such that xi € Bk for all I £ Ik- Therefore, 
in view of the fact that |ii| > 1, we must have |ifc+i| — |ifc| = for some 
1 < k < M — 1, implying that Bk+i\Bk H fi + = 0. Thus, there exists a cluster 
of droplets, whose indices are denoted by Ik, which are within 0(£~"( 1+fecr )) 
distance of xq and are separated from all other droplets by 0(e _Q ( 1 + cr + feo ")) 
distance. 

Let us show that this contradicts the minimality of E for small enough e. 
Indeed, let us displace the droplets in Bk to the new locations x\ — xi+\(xi— Xi), 
with I G Ik, which represents a dilation of Bk by a factor of 1 + A relative to 5.;, 
keeping all fi fixed. For < A <C 1 the resulting change AE of energy satisfies 

| In e| 2 A/5 < -cX + CXs aa \]iie\, (3.61) 

for some C, c > independent of e 1, where we used Lemmas |3.3|and|3.4[ and 



the estimate (2.4), arguing as in the derivation of (3.421. Thus, the considered 



rearrangement lowers the energy. □ 

As a simple corollary to this result, we actually have the following universal 
(^-independent) upper bound on \&f\ and, hence, on f^: 

Corollary 3.17. Let f2 + = LifL^tif be a non-trivial minimizer of E, where Ctf 
are the disjoint connected components ofCl + . Then, for any S > 



< 7r(l2(V2- l)) 2/3 + (5, (3.62) 



when e is sufficiently small. 

Proof. If is bigger, split Clf into two disks of equal area and move them 
apart a distance d = e - ' 3 , with < (3 < a < |. Arguing as before, the energy 
change AE upon this manipulation is given by 

| In el AE < 2(\/2 - 1)Vtt \n+\ 1/2 - ^\n+\ 2 + o(l). (3.63) 

ZTT 
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In view of the arbitrary closeness of (5 to |, the energy change is, therefore, 
negative for e -C 1. □ 



Finally, we note that the argument of Proposition |3. 16 still holds for local 
minimizers of low energy, the result can is obtained by sending A — > in the 
proof. 



3.3 Limiting behavior 

We now investigate the limiting behavior of the minimizers of E as e — > 0, with 
8 > ^\/9k 2 fixed, i.e., the situation in which minimizers are non-trivial. As 
the value of e is decreased, the number of droplets in a minimizer are expected 
to grow. What we will show below is that in the limit e — > the droplet 
sizes become asymptotically the same, and that the droplets become uniformly 
distributed throughout Q. 

Let us first study the behavior of each droplet as e — > 0. We have the 
following result. 



Proposition 3.18. Let f2 H 



flf be a non-trivial minimizer of E, where 



Q~f a re the disjoint connected components offl + , and let r, be as in Proposition 
Then fi — > y3 uniformly as e 



3.13 



0. 



3.14 



we already know that fj > \r3 — 5 
Let us prove that the matching upper 

p(xi) € Cl be 

a disk of radius centered at Xi defined in Pro positi on 3.13 and consider 
Vtp = H\ l JfL 1 B e -/s(xi). Note that by Propos ition 13.16 the disks B e -0(xi) do 



Proof. First of all, by Proposition 
for any S > 0, provided that e <C 1 

bound also holds for e <C 1. Indeed, for any j3 S (0, |) le t B F 



not intersect for e<l. In fact, by Proposition 3.16 for any a S |) we have 
dist{B e ~i3{x i ),B e -t,{x j )) > e~ a for e < 1. 

Let us show that the minimum of v defined in (3.4 1 is attained in tip for 
£ <C 1, Let x be such that v(x) = min and let Xi be the center of a droplet 
which is closest to x. Recalling the definition in (3.41 ) and Proposition 3.13 
can write 



v(x) 



Vi(x) 



2|lne 



ln( £ 1 / 3 (l + |5-^|))+o(l), 



(3.64) 



where we used (|2.4| 



In particular, for any 6 > we have v(x) > Vi(xi) + 
Xi\ < and e <C 1, in view of (3.42), where according 

On the 



3/3) - 5, if \x 

to Proposition |3.16[ we can use a defined above, whenever e <C 1 
other hand, choosing 7 € (/?, a) and picking any a;' such that \x' — Xi\ = e 
we see that for any S > we have w(x') < Vi(xi) + \r^(\ — 37) + 8 for e«l. 
However, with 5 sufficiently small this implies that v(x') < v{x) for small enough 
e, contradicting minimality of v at x. 

Now, we demonstrate that v(x) > \5kT 2 — |v9 — 5, for any 6 > 0, provided 
that e«l. Indeed, suppose the opposite inequality holds for some 6 > and a 
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sequence of e — > 0. Then, inserting a new droplet in the form of a disk of radius 
f = 0(1) centered at x results in the change AE of energy 



\lne\AE = V v( x)(r) + o(l), 



(3.65) 



where V is given by (3.181, and we used (2.4 1 and (3.42). Since by assumption 
2v(x)-Sk- 2 < |^9, it is easy to verify that V v /^ attains a minimum at some 
f = f o > s/3, with Vu( 2 )(fo) < 0. Therefore, inserting a droplet with radius fo 
and center at x would reduce energy for some e< 1, contradicting minimality 
of E. 

This, in turn, implies that Vi(xi) > \5n^ 2 — | \/9 — <5 for all i. Indeed, since 
x £ tip, from ( Egl we have > ifc" 2 -^^9- ^(1 - 3/3) + o(l), for any 

x' € dB e -p{xi), for E<Cl. On the other hand, by (3.42) the same inequality 
holds for x! — Xi. The estimate then follows in view of arbitrariness of j3 < |. 

Finally, since ^i(xj) > \8kT 2 — j\/9 — S when e -C 1, and by Proposition 
3.13 li) the values of fj are close to the minimizers of V v .^ s .^(f) for r > v3 — 5, 
by direct inspection we have f j < v^3 + (5 as well, for any 5 > and e <C 1. □ 



Let us point out that by (3.461 the uniform convergence of the droplet 



radii in Proposition 3.18 also implies uniform convergence of Vi(xi) to a space- 
independent constant as e — > 0: 




(3.66) 



In fact, from the pr oof o f Proposition |3.18| we can conclude that v stays close 
to the constant in (3.661 in tip for (3 arbitrarily close to |, provided e -C 1. 
This, in turn, implies that the droplets become uniformly distributed in f2, or, 
equivalently, in Q, as e — > 0. Below we prove this fact, which also gives the 
leading order behavior of energy in the limit. 

Let us rewrite the energy E^ for the system of interacting droplets, using 



(3.18): 



1 N 

En = 7j^| E (V Vz ( Sz )(n) - £nvi(xi)i 



»=i 



4^ 2 
liner 



iv-i w 



£ E G( £ 1 /3| ln£ |-l/3 ( ^_ 2 . ))f 2 f | +o(1) . 



(3.67) 



i=l j=i+l 



To proceed, let us go back to the original scaling in x and introduce Xi = 
e 1 / 3 ! In el -1 / 3 ^. Also, for any < a -C 1 define (our method is reminiscent of 
the Ewald summation technique) 



G a {x) 



47T 2 £ 2 



/ e ^ if (K|j/-n|)dj/. 



(3.68) 
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Here G a is a mollified version of G. with Fourier transform 



G a (q)= / e^ x G a (x)dx 



K 2 +\q\ 2 ' 
and which can, e.g., be estimated as 

G„{x) = G(x)+o(e< J ' 4 ), \x\>e a ' 2 , 

and 

G (T (x) = 0(a\lne\), \x\ < e a . 
Therefore, in view of Lemmas |3.4| and |3.6| we can write 

1 N 

En = [kid E ( Vv ^( f ^ ~ 47r Vi(xi)r?) 



(3.69) 



(3.70) 



(3.71) 



»=i 



2tt 2 
iTnel 



(3.72) 



i=l 3=1 

Now, let us introduce the quantity 



1 " 

p( a ;) = _- r 2 < J( a ;- a:< ) 



(3.73) 



Note that by L emma 3.6 we have J Q p(x)dx = 0(1). In view of Proposition 



3.18 and (3.661, we can further rewrite En as 

A? 



E 



N 



2ttv/9 



^9 



111).- ^ •' "-A r ^ K 2 I 5 2 



p{x)dx 



+6?r 2 v / 3 / / p{x)G a {x-y)p{y)dxdy + 0{<j). (3.74) 



In fact, by Proposition 3.18 and (3.66), the first term in (3.741 goes to zero. 

(3.75) 



Therefore, in terms of the Fourier coefficients 

Pq — i e lq ' x p(x) dx 
we can write 



En 



2ir^9 - s/9 



+6tt 2 ^3 



gS27rZ 2 \{0} 



K ] p + 



e^'\"\ 2 \p q \ 2 
n 2 + \q\ 2 



6tt 2 ^3 



0(a). 



(3.76) 
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Minimizing this with respect to po, we obtain 

2 



E N > 



1 

2^2 



^9 



6tt 2 ^3 

? e2irZ 2 \{0} 



e 2 



\h\ 2 



0(a). (3.77) 



Finally, from Lemma 



A.2 



one can see that min E N < — ^ (3 — | \/ 9 k 2 ) + 
o(l) for £ < 1. Hence, in view of arbitrariness of a the constant in (3.771 is 
the limit of Em and, by Proposition 3.15 also of E, as e — > 0. In addition, this 



implies that p q — > as e — > for every q ^ 0. Thus, we just proved 



Proposition 3.19. Let 6 > k 2 , and let p be defined in with Xi 

e 1 / 3 ! In e\~ 1 /^Xi, where Si are as in Proposition 



3.13 



Then 



min E 




and 



^9 



2vr^ 



(3.78) 



(3.79) 



weakly in the sense of measures, as e — > 0. 

We end by noting that the homogenization approach to multi-droplet pat- 
terns in a related context was first discussed in [57]. Also, let us mention that in 
a related class of problems existence of limiting density for the ground states of 
particle systems interacting via potentials like our G as the number of particles 
goes to infinity was proved in [58]. The difference with our result, however, is 
that in [58] the limit is taken at fixed positive temperature, while in our case the 
system's temperature (in the usual thermodynamic sense) is strictly zero. Yet, 
as was pointed out in [58], the "effective" temperature of the system considered 
actually goes to zero as the number of particles goes to infinity, making these 
results closely related to ours. 



3.4 Fine structure of the transition point 

Finally, we briefly discuss the appearance of non-trivial minimizers in the vicin- 
ity of the point S m = ^\/9k 2 (this transition point was identified in [13,14]). 
First, we note that by the results just obtained the transition from trivial to 
non-trivial minimizers appears to be quite abrupt in the limit e — > 0. In fact, 
in this limit one goes immediately from no droplets to infinitely many droplets 
upon crossing the point 8 — 5 m from below. 

Observe that the energy E[u) (and, equivalently, E[u] — E[— 1]) is amonoton- 
ically decreasing function of 5. Therefore, a passage through the neighborhood 
of 5 = S m at small but finite e will result in a monotonic increase of the number 
of droplets in a minimizer. This number will quickly get large as one moves 
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away from the transition point. Therefore, in order to analyze droplet creation 
at the transition, we need to further zoom in on the parameter region around 
S = 5 m . Let us introduce the renormalized distance to the transition (with the 
transition point shifted appropriately): 



I In el 




(3.80) 



and consider the behavior of energy E in the limit e — > for r = 0(1). As can 
be easily seen, all the estimates obtained previously remain valid in this case, 
and minimizers are close to a collection of N disks separated by large distances, 
whose energy is given by Em to 0(e a ). We can also write the energy Em in the 
form 

N N-l N 

^ = ^£ 1 (f i )+4 7 r 2 |ln£|- 2 ^ ^ G( Xi - Xj )r*r], (3.81) 

i—1 i—1 j— 2+1 



where 



|lne|Si(f) = 2tt Or - 



f 2 + if 4 



+ |7r| lnel" 1 In | lne| (f 2 - ^9)f 2 

— tt[ lnej- 1 ^ 2 ^^ — |) + 2r)r 2 (3.82) 

is the energy of one disk-shaped droplet of radius f. 

It is easy to see that in the limit e — > we have E\(f) > for all f > 0, 
and E\{f) = if and only of f = v3. Therefore, fj — ► s/3 uniformly in a 
minimizer as e — > with r fixed. In fact, by convexity of E\ near f = v3 we 
have fj — -^3 = 0(\ lne| _1 In | In e|) in the limit e — > 0. Therefore, we obtain 
(the summation is absent in the formula, if N = 1) 



E N = 12\/3n 



o(|lne|^ 2 ). (3.83) 

From this expression it is easy to see that N — O(l) quantity. Thus, in this 
case the problem reduces to minimizing the pair interaction potential given by 




the sum in (3.83). We summarize the above discussion by stating the following 
result. 



Proposition 3.20. Let S be given by (3.80) with r fixed. Then, there exists 



a strictly monotonically increasing sequence of numbers (r n ) , with r„ — > oo as 
n — > oo, such that, provided that s <C 1." 

(i) If t < t% = ^jv^9(3 — 41n3 — 121nS), then there are no non-trivial 
minimizers of E. 
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(ii) If Ti < t < T2, with Ti = T\ + 2n-\f§ min G, the minimizer of E is a single 
droplet. 

(Hi) If T n < t < T n+ \, all minimizers of E consist of precisely n droplets. The 

to— 1 n 

droplet centers {x;} nearly minimize V — G(xi — Xj). 

i=l j=i+l 

Let us mention that local minimizers of E without screening (i.e. with 
k — > 0) which are close to disks of the same radius centered at the minimizers 
of V were constructed perturbatively in a recent work of Ren and Wei [32,33]. 
We note that when t = 0(1), existence of these solutions easily follows from 
our analysis, if one notices that in the considered regime the excess energy of 
a minimizing sequence controls the isoperimetric deficit of each droplet and 
enforces 0(1) distance between them. Therefore, solutions with a prescribed 
number of droplets may be obtained by minimizing over all u G BV{£L\ {—1, 1}), 
such that the support of 1 + it has a fixed number of disjoint components. In 



turn, by Proposition 3.20 the global minimizers of E belong to this class. 



4 Connection to the diffuse interface energy 

We now turn to the study of the relationship between the sharp interface energy 
E and the diffuse interface energy £. Since most of our analysis here does not 
rely on any particular assumptions on the dimensionality of space, we will treat 
the general case of f2 being a d-dimensional torus: f2 = [0, l) d . We assume that 
W is a symmetric double- well potential with non-degenerate minima at u = ±1, 
together with some natural technical assumptions: 

(i) W G C 3 (R), W{u) = W(-u), and W > 0, 

(ii) W(+l) = W(-l) = and W"(+l) = W"{-1) > 0. 

(iii) W"(|u|) is monotonically increasing for |u| > 1, lim^^^W" (u) = +oo, 
and \W'(u)\ < C(l + \u\ q ), for some C > and q > 1, with q < if 
d>2. 

Since we are setting the surface tension to e, we need to additionally normalize 
W as follows: 

(iv) We have 

y/2W(u) du=l. (4.1) 



Note that these assumptions are satisfied for, e.g., the rescaled version of the 
classical Ginzburg-Landau energy: W(u) = J|(l — u 2 ) 2 for d < 3. Also note 



that this assumption is not restrictive, since it is always possible to make (4.1 1 
hold by an appropriate rescaling. 
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Let us begin our analysis with a few general observations. First of all, it is 
clear from standard arguments (see e.g. [59]) that for any e > there exists a 
minimizcr u £ H l (£l) of £ satisfying J^udx — u. Note that any critical point 
u of £, including minimizers, is a weak solution of the Euler-Lagrange equation 



(here and below Go solves (1.2 1 with periodic boundary conditions and has zero 
mean) 



e 2 Au - W'(u) - v + n = 0, v = I G Q (x - y)(u(y) - u)dy, (4.2) 

Jn 

where 

(i= ( W'{u)dx (4.3) 
Jn 

is the Lagrange multiplier. Furthermore, from the Sobolcv imbedding theorem 
we have u £ for p = and hence v £ W 2 ' p {Cl) C C°< a (f2), for some 

a £ (0, 1), if d < 6. Applying Moser iteration [59], we then find that u £ L p (£l), 
for any p < oo. Therefore, by standard elliptic regularity theory [50], we also 



have u £ W 2 ' p (£l), so u £ L°°(f2) and is, in fact, a classical solution of (4.2 1. 

We now show that u is uniformly bounded and that |u| cannot much exceed 
1 whenever £[u] is sufficiently small, at least for d not too high. 

Proposition 4.1. Let d < 6 and let u be a critical point of £. Then, for every 
5 > we have \u\ < 1 + 5 and \v\ < S in f2, whenever £[u] is sufficiently small. 

Proof. Observe first that for every 5 > and £[u] small enough we have |{|w| > 
1 + S}\ < |. Now, suppose that the maximum value u m of |u| is greater than 
1 + 5. Without a loss of generality, we may assume that u m = maxu. By the 
preceding observation, we have \i < C + \ W{u m ), for some C > independent 
of u rn . Therefore, in view of the monotonic increase to infinity of W'(u) due to 
hypothesis (iii) on W, we have [i < ^W'(u m ) for u m sufficiently large. Now, 
taking into account that v > —C'u m for some C > and large enough u m , from 



(4.2 1 we find that e 2 Au > \W{u m ) — C'u m > at the point where u = u n 
in view of assumption (iii) on W, contradicting the maximality of u. Finally, 
to see that |u| < 1 + 6 with any 8 > 0, when £[u] <C 1, note that u — > ±1 a.e. 
when £[u] — > 0. Hence, in view of the uniform bound on u obtained earlier, we 
have fj, — » 0. Furthermore, since the non-local term in the energy can be written 
as \ J n \Vv\ 2 dx, and v is uniformly bounded in W 2,P (Q) for any p < oo, we 
also have v — > uniformly in f2 in this limit. Therefore, in view of positivity of 
W'(l + S), we can apply the same argument as above to complete the proof of 
the statement. □ 

We note that while the arguments above hold for every critical point £ with 
small energy, it is generally possible for a local minimizer of £ to strongly deviate 
from ±1 in most of f2: take, for instance, one-dimensional periodic solutions of 



(4.2 1 with period O(l) [4,60,61]. Of course, these critical points will have O(l) 



energy when e — > 0, as opposed to minimizers of £ whose energy vanishes in 
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this limit. Let us also mention that numerical evidence shows that generally 
max | it | > 1, even for minimizers and e <C 1. 

We now turn to estimating the minimal energy of £ from below by the 
minimal energy of E. For u G with J udx — u, let us separate the 

domain f2 into three pairwise-disjoint subdomains: 

fi = U fti U fig. (4.4) 

where 

Q 5 + = {xeQ: u(x) >l-8}, (4.5) 

n s _ = {xen-.u(x)<-l + 5}, (4.6) 

Q 5 = {xeQ:-l + S<u{x)<l-6}, (4.7) 

Next, let us introduce the following three auxiliary functionals (for simplicity of 
notation, we will suppress the index S in the definition of each functional): 

Si[u] = J^(^\Vuf + W{u))dx, (4.8) 
£ 2 [u] = / (u-u a ) 2 dx 



2k 



Q S + UQ S _ 



+ If [ \u{x)-u)G {x-y){u{y)-u)dxdy, (4.9) 
1 Jn Jn 



V5 



where uq(x) = ±1 whenever x € f2j_, respectively, with 

1 



and 

1 



(4.10) 



[«]= / 



W(u)-^2 («-«o) 2 )<to- (4.11) 



It is clear that the energy £ can be estimated from below as 

£ > £i + £i + £3. (4.12) 

Hence, we are going to establish a lower bound for £ by considering the lower 
bounds for each term in the sum above. 

We start with the part of energy that is associated with the interfaces: 

Lemma 4.2. Let 5 > be sufficiently small, let u e H 1 ^) and suppose that 
\Q S _ U fl s + \ > 0. Then there exists u £ BV(Q;{-1,1}) such that u (x) = ±1 
whenever x G Q±, and 

£i[u] > |(1 -ai<5 2 ) / \X7u \dx (4.13) 
2 Jn 

for some a\ > independent of 5 and e. 



32 



Proof. First of all, if either \Q s + \ or \Q, 5 _\ is zero, we can simply choose uq to be 
constant (e.g. it = —1 when |0^_| = 0). So, let us assume that both \Cl±\ > 0, 
and approximate u in H l {£l) by a piecewise linear function u with Vu ^ 
almost everywhere in fL Then, using the Modica-Mortola trick [62,63] and the 
co-area formula [64] , we find 



£i[u]>e [ y/2W(u) |V«| dx = e [ y/2W(t) \{u = t}\ dt. (4.14) 

J\u\<l-6 J-l+8 



Since the function — t}\ is continuous for all t 6 [— 1 + 5, 1 — S), there exists 
a constant c e [— 1 + 5, 1 — <5] such that the right-hand side of (4.141 equals 

| {it = c}| f2~* d V 2W (*) dt > (! - M 2 )|{u = c}], for some a x > and all 5 
small enough. 

Now, define u G W(fi; {-1, 1}) as 

~ ( \ J +1 ' > C ' ^ Kl 

I — 1, u(x) < c. 

The preceding arguments imply the desired inequality for u. Passing to the 
limit in the approximation, we obtain the result, with uq — limuo in L x (0) 
upon extraction of a subsequence. □ 



Lemma 4.3. Let u and uq be as in Lemma let u satisfy J^udx — u, and 
let \u\ < 1 + 5 3 and \J n Gq(x — y)(u(y) — u) dy\ < <5 3 in fl, for 5 > sufficiently 
small. Then 

£%\u] > \ [ [ (u (x) - u)G(x - y)(u (y) - u) dxdy - a 2 6 £[u], (4.16) 
* Jn Jn 

for some a 2 > independent of 5 and e, whenever £[u] < S^ d+6 \ 
Proof. Let us write u as follows 

u = uo + ui + U2, ui(x) = -k 2 Go(x — y)(u(y) — u)dy. (4.17) 

Jn 

Note that by assumption, || William) < CS 3 and 1 1 Ua 1 1 z.°° (£1) — f° r some 
C > 0. Now, observe that u\ solves 

— Alii + K 2 Ul = —K 2 (uq + 1*2 — u), (4-18) 

and, therefore, we also have 

Ui (x) = -k 2 G(x-y)(u (y)+u 2 (y)-u)dy. (4.19) 
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Substituting u in the form (4.171 into (4.9), we obtain 
S 2 (u + ux + u 2 ) = tt^ [ (u 1 +u 2 ) 2 dx 



n*_un s _ 



(u (x) + ui(x) + u 2 (x) - u)G(x - y)(u Q (y) + u 2 (y) - u)dydx 



inJn 

^Ls uidx -hL uiu2dx -l 



u 2 (x)G(x - y)u 2 (y)dydx 



n Jn 



l, [ I ( u o(x) - u)G(x - y)(u (y) - u)dydx + 
z Jn Jn ZK Jn 5 + un s _ 

1 



u 2 dx 



- \ i f ( u o0) - u)G{x - y)(u (y) - u)dydx - I u{ dx 
2 Jn Jn 2k 2 J n s 

+ / / u 2 (x)G(x - y)(u (y) - u) dydx 
Jn s Jn 

+i / I uj(x) - k 2 u 2 (x) I G(x - y)u 2 (y)dy) dx. (4.20) 

2k Jn s + unt \ Jn s + un s _ J 



In fact, the last line in (4.201 is non-negative. Indeed, writing the integral in the 



last line of (4.20) with the help of the Fourier Transform d q of it = u 2 x^s un s , 

•is 1 1 r><5 . 



where Xn s + un s 1S the characteristic function of fl + U ft 

d q = e lq ' x u 2 (x) dx, 

Jn 6 + un s _ 

we obtain 



(4.21) 



u\{x) — k 2 u 2 (x) 



n s + un 6 _ 



G{x - y)u 2 (y)dy^j dx 

- E ™>o. 

q£ 2 TTl 



k 2 + |gf 



(4.22) 



To estimate the remaining terms in (4.201, we note that 

~2 [ Uiu 2 dx+ I I u 2 (x)G(x - y)(u (y) - u) dydx 
K Jn* Jni Jn 



In* Jn 



< 



nf, J {\u\>i-83} 



n s J{\u\<i-63} 



u 2 (x)G(x - y)u 2 {y) dydx 
\u 2 (x)\G(x - y)\u 2 {y)\dydx 
\u 2 (x)\G{x-y)\u 2 {y)\dydx. (4.23) 
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Since G £ L P (Q) for all p < (any p < oo in d = 2), by Holder inequality we 
can see that for any Cl C il 

G(x-y)h 2 (2/)|dy<C^K|Mx) ^ < C| |« 2 ] (4.34) 



for any ? > f • Therefore, continuing the estimates in (4.23), we obtain 



ni Jo 



u 2 (a;)G(a: - y)u 2 (y) dydx 

< c (V + s- 6/q £ 1/q [u}^ \n s \ < 2cs 3 \n s \, 



(4.25) 



whenever £[u] < 8 3 ^ q+2 \ where we took into account that by the assumptions 
of the lemma \u 2 \ < \u — u \ + \u%\ < CS 3 in {|it| > 1 — <5 3 }, and that £[u] > 
/{|«|<i-<5 3 } w ( u ) dx ^ «5 6 |{M < 1 - S 3 }\ for some c> 0. 
Similarly, we have 



/ u\ dx 


+ 


/ U1U2 dx 


< cs 3 \n 5 \. 


(4.26) 













The statement of the lemma then follows from the fact that £[u] > J n s W(u)dx > 
c<5 2 1 f^Q I , for some c > 0, by choosing q = |(d + 2). □ 

Lemma 4.4. Let u and ug be as in Lemma \4-2\ Then 

£ 3 [u] > -a 3 5£[u]. (4.27) 
for some a 3 > independent of 5 and e, for sufficiently small 6 > 0. 

Proof. By assumption (iii) on W, we have W(u) > 2^2 (u — uq) 2 whenever 
\u\ > 1. Hence 

1 



(4.28) 

□ 



W(u)- —(u-u o y\dx 



£a[«] ' / , , , 

J {1-S<\u\<l} V Zk 

(u — uo) 2 dx > — a 3 5£[u 



> -CS 



{1 — <5<|u|<l> 



for some 03 > 0. 



Combining all the results above with an observation that \Q S + U Q s _ \ > for 
£[u] small enough, we obtain 

Proposition 4.5. Let S > be sufficiently small, letu S H 1 ^) satisfy J n udx — 
U, let \u\ < 1 + S 3 and \J^Gq(x — y)(u(y) — u) dy\ < S 3 in ft, and let £ < 
£§(<i+6)_ Then there exists a fu nction uo E BV(Cl; {— 1, 1}) such that £[u] > 
(1 — 5 1 ' 2 )E[uq], with k given by (4-10). 
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Importantly, the lower bound in Proposition |4.5| is sharp in the limit e — ► 
for all functions uo G BV(Cl; {— 1, 1}) obeying suitable bounds (satisfied by 
minimizers of -E in d = 2) : 

Proposition 4.6. Le£ u € SV^fi; {— 1, 1}), iwift ifte jump se£ o/ class C 2 , 
let the principal curvatures of the jump set of uq be bounded by e~ a for some 
a € [0, 1), let the distance between different connected portions of the jump set 
be bounded by e a , and let |/ a G(x — y)(uo(y) — u) dy\ < S for some S > small 
enough. Then there exists a function u G iJ x (f2) with J n udx — u, such that 
< (1 + ^ 1 ^ 2 )E[uo], with k given by (4-10), whenever E[uq] < <5s( d + 3 ) and 
e < 1. 

Proof. For simplicity of presentation, we only give the proof in the case d = 2. 
With minor modifications, the proof remains valid for all d. 

Here we adapt the standard construction of a trial function for the local 
part of the Ginzburg-Landau energy. Let U(p) be the solution of the ordinary 
differential equation 



d 2 U 
dp 2 



W'(U) = 0, U(-oo) = 1, U(+oo) = -1, U(0) = 0, (4.29) 



where the last condition fixes translations. As is well-known (see e.g. [65]), 
this solution exists, is unique and is a strictly monotonically decreasing odd 
function, approaching the equilibria at p = ±oo exponentially fast. Therefore, 
for any 5 > we have \U(p)\ < 1 — 8, it and only if \p\ < /, with some positive 
/ = 0(| In <5 1 ) . Also note that by hypothesis (iv) on W 




dU 



dp 



W{U) } dp - 



1-5 



^2W(s) ds= l + 0(5 2 ). (4.30) 



1+4 



Now, introduce the signed distance function r(x) — ±dist(a;, CI ), where 
fl ± = {uq = ±1}, whenever x € CP, and define a regularized version Uq of uq: 



'U{e~ l r{x)), 



\r(x)\<el, 



u e (x) 



(l-S + e-^drix^-el^uoix), el < \r(x)\ < e{l + 1), (4.31) 

U {x), 

Then, it is easy to see that the function 



\r(x)\>e(l + l). 



u(x) = u%{x) - k 2 / G(x - y){u e (y) - u)dy. (4.32) 
Jo. 

is in H 1 ^), with J^udx — u. Moreover, we have for any q > 1 

\u(x) - ul(x)\ < k 2 / G{x - y)(u (y) - u) dy 
Jn 

G(x - y)\u £ (y) -u (y)\dy 
< C(6+ ICltl 1 ^) <C'{5+\\n8\ 1/q E 1/q [u ]) <C"S, (4.33) 



3G 



where we defined fli — {\r\ < e(l+ 1)}, estimated f n G{x — t/)(tig(y) — Uo(y)) dy 
as in Lemma 4.3 and used the curvature bound and the assumption on E with 
d = 2 and q = 2. On the other hand, after a few integrations by parts, from 
(4.321 we also obtain 

|V(it — Uq)\ 2 dx — — k 2 / (u — Uq)(u — u) dx 
Jn 

(u(x) — u)Gq(x — y)(u(y) — u) dydx < 2k 4 £[u]. 



K 4 



(4.34) 



To estimate £[u], let us introduce a system of curvilinear coordinates (p, £) 
consisting of the signed distance p to the jump set of Uq and the projection £ 
onto the jump set. By our assumptions this is possible whenever |r(x)| < e a . 
Therefore, for £<lwe can write 



£[u] 







du 




I (t 

Jon+ \ z 


dp 



du 



W(u) 



x(l + pK) dH 1 (£)dp+ [ ( e —\\7u\ 2 + W{u)\dx 
Jn\n, \ 2 J 



(u(x) — u)Gq(x — y){u{y) — u) dydx, (4.35) 



n Jn 



where K = K(£) is the curvature at point £ on the jump set of uq. Substituting 
the ansatz of (4.321 into (4.35), taking into account that |V(w — u§)| < C in 
fli for some C > independent of e (we have u — Uq uniformly bounded in 
W 2,P (Q), for any p < oo) and that V« = V(u — uq) = V(u — u§) in ri\fi;, and 
using (4.331 and (4.341, we obtain (estimating each line in (4.351 separately) 

£[u] = ^(l + 0(E 1 - a \la5\) + 0{d\lnS\)) [ \Vu \dx 
2 Jn 

(u - u ) 2 dx + 0(e 2 £[u}) + 0(6 £[u}) 



1 

~2^2 
1 

+ 2 



a\nj 



(u(x) - u)G (x - y){u(y) - u) dydx. 



(4.36) 



n Jn 



Now, using the identity 

(u(x) — u)Go(x — y)(u(y) — u) dydx + n~ 2 I (u — u £ Ydx 



n Jn 



(uq(x) - u)G{x - y){ul(y) - u) dydx, (4.37) 



n Jn 
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we can further write (4.361 as 



E[u] = E[u ] + 0(5\ln5\E[uo\)+0(5E[u])--^ / (u-u £ ) 2 dx 



ill 



(u 6 (x) - u (x))G(x - y)(u E (y) - u) dydx 
(u e (x) - u (x))G(x - y)K(y) - uo(y)) dydx, (4.38) 



n, Jn 



for £«1. Finally, using the same estimates as in (4.331, we obtain 

(1 + 0(5))£[u] = (1 + 0(S\ ]nS\))E[u ] + 0(6\Sl t \) + 0(|^| 3/2 ) 

= (l + O(5\]nS\))E[u ], (4.39) 

from which the result follows immediately. □ 

The last two propositions show asymptotic equivalence of the diffuse inter- 
face energy £ with the sharp interface energy E for sufficiently well-behaved 
critical points and e < 1. In particular, the energies of minimizers of both E 
and £ are asymptotically the same in the limit e — > 0. It would also be natural 
to think that the minimizers (even local, with low energy) of £ are, in some 
sense, close to minimizers of E when e <C 1 (this will be a subject of future 
study) . 



5 Proof of the theorems 

Here we complete the proofs of Theorems |2.1f|2.3| 

Proof of Theorem |2.1| The main point of the proof is the lower bound in 



(2.1 1, since the upper bound is easily obtained by constructing a suitable trial 
function (as in Lemma A.l I. The basic tool for the lower bound is a kind of 
interpolation inequality obtained in Lemma B.l Note that the proof for E 
works in any space dimension. 

To prove the lower bound, let us denote by u a minimizer of E. Introducing 



e lq - x {u{x) -u)dx, (5.1) 

Jn 

where q £ 2irZ d , we can estimate the energy of the minimizer as follows 



mini? > - [ I (u(x) — u)G(x — y)(u(y) — u) dxdy 
2 Jo Jn 

ij2 > M 2 



2 ^ K 2 + \q\ 2 ~ 2k 2 



2 k 



X -A f (u-u)dx\ = ^ 2 (\n+\- l -±^\ , (5.2) 
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where we introduced the set f2 + = \u = +1} 



In view of the upper bound in (2.1 1, it follows from (5.2 1 that |f2 + | = |(1 + 
u) + 0(e 2 / 3 ), implying that is bounded away from or 1 for e <C 1. Hence, 
by isoperimetric inequality there exists p > such that 



P 



|Vl(| > p, 



(5.3) 



whenever e< 1. Applying Lemma |B.1| to u — u, we conclude that 

C 

vamE > eP + ^ , 

for some C > independent of e, for e -C 1. The result then follows from an 
application of Young inequality and Propositions |4.1| and |4.5| □ 



(5.4) 



Proof of Theorem |2.2| This theorem combines a number of results proved 
in Sec. [3] in the original, unsealed variables. Part (i) of the theorem is the 
statement of Prop osition 3.5 Part (ii) of the theorem is the collection of results 

fc ] < E[-l] = l £ 4 /3| ln£ |2/3 K -2^2 

and 



from Lemma A.l (taking into account that E[v,k 



for 5 > 



3.16 



Corollary 



3.2 



Lemma 



3.6 



and Propositions 



3.13 



3.15 



with a = i — a. Part (hi) of the theorem is contained in the statements of 



Propositions |3.18| and |3.19| 



□ 



Proof of Theorem 2.3 



First of all, we have min£ <ti 1 when e <C 1 and 
since min£ < £{u) = 0(e 4//3 | In e| 2 / 3 ) in that case. 



u = -1 + 0(e 2 / 3 |lne| 
Then, from Proposition |4.1| and Lemma |3.8| we conclude that the assumptions 
of Propositions |4.5| and |4.6| are satisfied for the minimizers of £ . Therefore, 
the energies E and £ are asymptotic ally t he same in the considered limit, and 
the conclusion follows from Theorem 
monotone decrease of E with 6) . 



2.2 



(the case 5 — k \/9 n 2 is included by 



□ 
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A Upper bound 

Here we construct a trial function that achieves the lower bound for the energy 
of the non-trivial minimizers of E. 

Lemma A.l. Let u = — 1 + e 2/>3 \ lne^^S, with 5 > \\/%k 2 fixed. Then there 
exists u G BV(fl; {—1, 1}), such that 

E[ „] _ £ </a| her { f (J - \ W>« 2 ) + O (^) } , (A.l) 
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for £«1. 

Proof. First, consider u\{x) = — 1+2xb t {o){ x ), where Xs r (o) is the characteristic 
function of a disk of radius r centered at the origin. If v\(x) — J n G(x — 
y)(u\(y) — u) dy, then by using (2.3 1 we explicitly have (see (3.40 1 ) 

M*) = + il 1 ~ KrKtWloiKlxfi), (A.2) 



2 



rh(Kr)K ( K \x + n\)), \x\<r, (A.3) 

K nGZ 2 \{0} 

where K n and /„ are the modified Bessel functions of the first and second kind. 
Therefore, expanding the Bessel functions for r < 1 [45], we can write for \x\ < r 

Ul ( a; ) = _i + ^_^(21n/tr + 27-ln4-l)-^- 
K 2 2 

+r 2 J2 ^o(«l^ + n|) + 0(r 4 |lnr|), (A.4) 

neZ 2 \{0} 

where 7 ~ 0.5772 is the Euler's constant. Substituting this expression into the 
definition of E, after integration we get 

E[ui] = 27TCT- + i(l + u) 2 kT 2 - 2tt(1 + u)«T 2 r 2 

— 7rr (In Kr + 7 — ln2-i) + 7rr 4 ^ ^ (/c|n|) + 0(r 6 | lnr|). (A.5) 

n£Z 2 \{0} 

Now, consider a new test function 

k k 

U k (x) = -1 + 2 XB r (e 1 (fc 1 -i)+e 2 (fc 2 -i))(^) 5 ( A -6) 

k]=l fc 2 =l 

consisting of fc 2 disks of radius r arranged periodically in £1 (here ei and &2 are 
the unit vectors along the coordinate axes). We have 

E[u k ] = |(1 + u) 2 «T 2 + tt/c 2 (2er - 2(1 + u) K " 2 r 2 

-r 4 (ln K r + 7-ln2- |) + r 4 ^ ^o(k/c~ Vl))+0(fc 2 r 6 | lnr|). (A.7) 

nez 2 \{o} 



Approximating the sum in (A.7 1 by an integral: 

k- 2 



V K (nk- 1 \n\)= [ K a {K\x\)dx + 0{k- 2 \nk) 



i£Z 2 \{0} 



= 27r^ 2 + 0(fc- 2 lnfc), (A.8) 
and expanding for r< 1, we can further write 

E[u k ] = |(1 + u) 2 k~ 2 + nk 2 (2er - 2(1 + u)KT 2 r 2 

-r 4 lnr + 27TK- 2 r 4 k 2 ^ + Q(k 2 r 4 In k). (A.9) 
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We now substitute r — e 1 ^ z \ lne| 1 ' 3 \^3 into the expression above. Using also 
the definition in (2.7 1, we can write 



E[u k ] = e 4/3 | In e\ 2 / z (\n- 2 ~5 2 - 2tt^9| In e\~ X «T 2 (s - \ ^9 k 2 ) k 2 

+67r 2 ^3^ 2 |lner 2 fc 4N ) +0(e 4 / 3 |ln£|- 4 / 3 fc 2 lnfc). (A.IO) 



Finally, setting 



Ijngj 
2tt^9 



0(1), 



we obtain (A.l) with u = it/.. 



(A.ll) 



□ 



Let us also quote without proof a similar result concerning the upper bound 
for the reduced energy E^. 



Lemma A.2. Let u = -1 + e 2/3 \ lne\ 1/3 6, with 6 > \ \/9 k 2 fixed. The 



min En < — 



2 



O 



£ 4 / 3 ln|lne|' 
llnep/s 



(A.12) 



B Interpolation inequality 

Here we present the lemma that connects the non-local part of the energy with 
the interfacial energy via a kind of an interpolation inequality between BV(Q,), 
iJ~ 4 (f2) and L°°(f2), for functions bounded away from zero. 

Lemma B.l. Let u G BV(Q), where f2 = [0,l) d is a torus, and assume that 
m < |u| < M in fl for some M > m > 0. Let also |Vu| dx > p > 0, and 
let G solve \1.5^ in fl with periodic boundary conditions. Then there exists a 
constant C — C(d, n/p : m 2 /M) > such that 

[ [ u(x)G(x - y)u(y) dx dy > C ( [ \Vu\dx) . (B.l) 
Jn Jn \Jn J 



Proof. First, extend u periodically to the whole of W 1 . Then, introducing 
Xs( x ) — S~~ d \Bi\ ~ 1 x(S~ 1 x), where \ is the characteristic function of the unit 
ball B\ centered at the origin, we have 

/ / u{x)xs{x ~ y)u(y) dydx = f f u(x)u(x + Sy) dy dx 

Jn JR" I-diI Jn JBt 



> m 2 - 



MS 

Jb 



- [ [ [ \Wu(x + 6ty)\dtdydx>m 2 -Ms[\Wu\dx, (B.2) 
il Jn Jb 1 Jo Jn 
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where the inequality is obtained by approximating u by C 1 functions and passing 
to the limit. Therefore, choosing 

we obtain 

m 2 ^ 



n 



2 



/ / u(x)xs(x~y)u(y)dxdy = VxiMI^ 2 . ( B -4) 



where we introduced Fourier transform u g of w: 



u q = / e^ufc) tic, (B.5) 
Jn 

with q 6 2ir1 d . The Fourier transform x<5 of xs is, in turn, explicitly given by 

Uq) = (^) d/2r (^ + 1 ) J «/*W q V> (B ' 6) 

where Jd/2(%) is the Bessel function of the first kind and T(x) is the gamma- 
function. Now, applying Cauchy-Schwarz inequality, we obtain 



4 

m 



"„ 2 



<su P {^(g)( K 2 + M 2 )}EKI 2 
q « 

u{x)G{x — y)u{y) dx dy . (B-7) 



/n jq. 

Taking into account that J2 q \™q\ 2 — IMIi»(n) < M 2 an d that [45] 

yC2\K m M l p + \q\ o z )(\q\d) a , \q\S > 1, 
for some Ci^ > depending only on <i, we conclude that 

CS 2 < [ [ u(x)G(x-y)u(y)dxdy. (B.9) 

for some C > depending only on d, n/p, and m 2 /M. The result then follows 
immediately from (B.3l. □ 

Let us also make some remarks regarding a few extensions of these argu- 
ments. First, the same estimate holds true in the case where G is the Green's 
function of the Laplacian in f2 and u has zero mean. Note that in this case 
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the constant C in (B.l I becomes independent on p. The proof easily follows by 
passing to the limit k — > in the lemma. Another observation is that, actually 
for the considered class of functions a stronger interpolation inequality involving 
negative Sobolev norms holds. We give only the statement of the result, the 
proof follows easily by modifying a few steps in the arguments above 

Proposition B.2. Let u be as in Lemma \B.l\ Then 

[ u(l - A)"^ udx > C ( [ \Vu\dx] , (B.IO) 
Jn \Jn J 

for some C — C(d, p, m, M) > 0. 



C First and second variation 



Here we present the derivation of the first and second variation of E in d = 2, 
adapted from [14]. 

Lemma C.l. Let Cl + C f2 be a set with boundary of class C 2 and v be given by 
(3.4-)- Then, the functional E is twice continuously Gdteaux-differentiable with 



respect to C -perturbations ofdil + . Furthermore, the first and second Gateaux 
derivatives of E are given by f£gU and (3.7). 



Proof. Let a > 0, let p E C 1 (dCl + ), and let f2+ be the set obtained from f2 + by 
transporting each point of <9f2 + by ap in the direction of the outward normal. 
Note that for sufficiently small a the set <9f2+ is of class G\ in view of regularity 
of dCl+. Then, if E a = E(Cl+) and E = E(Cl+), from we have explicitly 



\lne\(E a -E) 



ap(x 



+ aK(x)p(x)) 2 + a 2 \Wp(x)\ 2 - 1 dH^x) 



(4v(x + rv(x)) - 28k- 2 ){1 + K(x)r) drdH 1 (x) 

on+ Jo 

r r rap{x) rap(y) 

+2|ln £ r 1 / / / / (l + K(x)r)(l + K(y)r') 
Jan+ Jon+ Jo Jo 

xG(e 1/3 | lne\- 1/3 (x + rv(x) -y- r'v{y))) dr'dr dH l (y)dH 1 (x) , 



(C.l) 

where K{x) is curvature, v(x) is the outward unit normal at x € d£l + , and we 
rewrote the integrals in terms of the curvilinear coordinates consisting of the 
projection cc of a point x € Ct to dti + and signed distance r = v{x) ■ (x — x), 
which is possible for sufficiently small a. Now, Taylor-expanding the integrands 
in the powers of r and integrating over r and r', after some tedious algebra we 
obtain that for any a S (0,1) it holds 



E„ 



E + a 



dE a 
da 



a=0 



fl2 cPEa 
~2~ "doY 



+ 0{a 2+a ) 



(C.2) 



a=0 



43 



where the derivatives are given by (3.6) and (3.7). In estimating the remainder 
term in (C.2 1 we took into account that v € C i,Q (f2) and the following estimate 
of the terms involving the convolution integral: 

r a p(y) 



on+ Jo 



< C 

< C 



ap(y) 



dn+ Jo 
a 2 - 



(G(£ 1/3 | lne\- 1/3 (x + v{x)r -y- v{y)r')) 

-G(e 1/3 | \ne\- x ' 3 {x - y))\ dr'dn\y) 
\x — y + v(x)r — v(y)r'\ 



In 



d(l+n\x-y\>Ma Jo 

< Ca 2 



dr'dU 1 ^) 

\x-y\ 

x-y\ 



dH\y) 



< C'a 2 \lna\, 



dQ+n\x-y\>Ma \ x V\ 



(C.3) 



for a<Cl, where M > is sufficiently large, and we used the series expansion 
of G [45]. 

Finally for every sufficiently small G 1 -perturbation <9f2+ of d$l + the distance 
from a point x € 9f2 + to <9f2+ is a C 1 -function, hence the formulas obtained 
above apply to all such perturbations. □ 



References 

[1] Bray, A. J.: Theory of phase-ordering kinetics. Adv. Phys. 43 (1994) 357- 
459 

[2] Landau, L.D., Lifshits, E.M.: Course of Theoretical Physics. Volume 8. 
Pergamon Press, London (1984) 

[3] Grosberg, A.Y., Khokhlov, A.R.: Statistical Physics of Macromolecules. 
AIP Press, New York (1994) 

[4] Kerner, B.S., Osipov, V.V.: Autosolitons. Kluwer, Dordrecht (1994) 

[5] Vedmedenko, E.Y.: Competing Interactions and Pattern Formation in 
Nanoworld. Wiley, Weinheim, Germany (2007) 

[6] Muthukumar, M., Ober, C.K., Thomas, E.L.: Competing interactions and 
levels of ordering in self-organizing polymeric materials. Science 277 (1997) 
1225-1232 

[7] DeSimone, A., Kohn, R.V., Muller, S., Otto, F.: Magnetic 
microstructures — a paradigm of multiscale problems. In: ICIAM 99 (Ed- 
inburgh). Oxford Univ. Press (2000) 175-190 



44 



[8] Choksi, R., Conti, S., Kohn, R.V., Otto, F.: Ground state energy scaling 
laws during the onset and destruction of the intermediate state in a type I 
superconductor. Comm. Pure Appl. Math. 61 (2008) 595-626 

[9] Choksi, R.: Scaling laws in microphase separation of diblock copolymers. 
J. Nonlinear Sci. 11 (2001) 223-236 

[10] Seul, M., Andclman, D.: Domain shapes and patterns: the phenomenology 
of modulated phases. Science 267 (1995) 476-483 

[11] Yu, B., Sun, P., Chen, T., Jin, Q., Ding, D., Li, B., Shi, A.C.: Self- 
assembled morphologies of diblock copolymers confined in nanochanncls: 
Effects of confinement geometry. J. Chcm. Phys. 126 (2007) 204903 pp. 
1-5 

[12] Kohn, R.V.: Energy-driven pattern formation. In: International Congress 
of Mathematicians. Vol. I. Eur. Math. Soc, Zurich (2007) 359-383 

[13] Muratov, C.B.: Theory of domain patterns in systems with long-range 
interactions of Coulombic type. Ph. D. Thesis, Boston University (1998) 

[14] Muratov, C.B.: Theory of domain patterns in systems with long-range 
interactions of Coulomb type. Phys. Rev. E 66 (2002) 066108 pp. 1-25 

[15] Care, CM., March, N.H.: Electron crystallization. Adv. Phys. 24 (1975) 
101-116 

[16] Emery, V.J., Kivclson, S.A.: Frustrated electronic phase-separation and 
high-temperature superconductors. Physica C 209 (1993) 597-621 

[17] Chen, L.Q., Khachaturyan, A.G.: Dynamics of simultaneous ordering and 
phase separation and effect of long-range Coulomb interactions. Phys. Rev. 
Lett. 70 (1993) 1477-1480 

[18] Nyrkova, LA., Khokhlov, A.R., Doi, M.: Microdomain structures in poly- 
electrolyte systems: calculation of the phase diagrams by direct minimiza- 
tion of the free energy. Macromolecules 27 (1994) 4220-4230 

[19] Ohta, T., Kawasaki, K.: Equilibrium morphologies of block copolymer 
melts. Macromolecules 19 (1986) 2621-2632 

[20] Bates, F.S., Fredrickson, G.H.: Block copolymers - designer soft materials. 
Physics Today 52 (1999) 32-38 

[21] Matsen, M.W.: The standard Gaussian model for block copolymer melts. 
J. Phys.: Condons. Matter 14 (2002) R21-R47 

[22] de Gennes, P.G.: Effect of cross-links on a mixture of polymers. J. dc 
Physique - Lett. 40 (1979) 69-72 



45 



[23] Stillinger, F.H.: Variational model for micelle structure. J. Chem. Phys. 
78 (1983) 4654-4661 

[24] Ohta, T., Ito, A., Tetsuka, A.: Self-organization in an excitable reaction- 
diffusion system: synchronization of oscillating domains in one dimension. 
Phys. Rev. A 42 (1990) 3225-3232 

[25] Glotzer, S., Di Marzio, E.A., Muthukumar, M.: Reaction-controlled mor- 
phology of phase-separating mixtures. Phys. Rev. Lett. 74 (1995) 2034- 
2037 

[26] Matsen, M.W., Bates, F.S.: Unifying weak- and strong-segregation block 
copolymer theories. Macromolecules 29 (1996) 1091-1098 

[27] Choksi, R., Ren, X.: On the derivation of a density functional theory for 
microphase separation of diblock copolymers. J. Statist. Phys. 113 (2003) 
151-176 

[28] Muratov, C.B., Novaga, M., Orlandi, G., Garci'a-Cervera, C.J.: Geometric 
strong segregation theory for compositionally asymmetric diblock copoly- 
mer melts. In: Singularities in nonlinear evolution phenomena and appli- 
cations. CRM Series. Birkhauser (2009) (to appear). 

[29] Miiller, S.: Singular perturbations as a selection criterion for periodic 
minimizing sequences. Calc. Var. Part. Dif. 1 (1993) 169-204 

[30] Ren, X.F., Wei, J.C.: On the multiplicity of solutions of two nonlocal 
variational problems. SIAM J. Math. Anal. 31 (2000) 909-924 

[31] Ren, X., Wei, J.: On energy minimizers of the diblock copolymer problem. 
Interfaces Free Bound. 5 (2003) 193-238 

[32] Ren, X., Wei, J.: Single droplet pattern in the cylindrical phase of diblock 
copolymer morphology. J. Nonlinear Sci. 17 (2007) 471-503 

[33] Ren, X., Wei, J.: Many droplet pattern in the cylindrical phase of diblock 
copolymer morphology. Rev. Math. Phys. 19 (2007) 879-921 

[34] Ren, X., Wei, J.: Droplet solutions in the diblock copolymer problem with 
skewed monomer composition. Calc. Var. Partial Differential Equations 25 
(2006) 333-359 

[35] Roger, M., Tonegawa, Y.: Convergence of phase- field approximations to 
the Gibbs-Thomson law. Calc. Var. PDE 32 (2008) 111-136 

[36] Albcrti, C, Choksi, R., Otto, F.: Uniform energy distribution for an 
isoperimetric problem with long-range interactions. J. Amer. Math. Soc. 
22 (2009) 569-605 

[37] Nishiura, Y., Ohnishi, I.: Some mathematical aspects of the micro-phase 
separation in diblock copolymers. Physica D 84 (1995) 31-39 



46 



[38] Choksi, R., Pclcticr, M.A., Williams, J.F.: On the phase diagram for 
microphase separation of diblock copolymers: an approach via a nonlocal 
Cahn-Hilliard functional. SIAM J. Appl. Math. 69 (2008) 1712-1738 

[39] Pctrich, D.M., Goldstein, R.E.: Nonlocal contour dynamics model for 
chemical front motion. Phys. Rev. Lett. 72 (1994) 1120 1123 

[40] Goldstein, R.E., Muraki, D.J., Petrich, D.M.: Interface proliferation and 
the growth of labyrinths in a reaction-diffusion system. Phys. Rev. E 53 
(1996) 3933-3957 

[41] Yip, N.K.: Structure of stable solutions of a one-dimensional variational 
problem. ESAIM Control Optim. Calc. Var. 12 (2006) 721-751 

[42] Ertl, G.: Reactions at surfaces: From atoms to complex- 

ity http: / / nobelprize.org/ nobeLprizes / chemistry /laureates /2007/ ertl- 
lecture.html (2007) 

[43] Ertl, G.: Reactions at surfaces: From atoms to complexity (Nobel lecture). 
Angcw. Chcm. Int. Ed. 47 (2008) 3524-3535 

[44] Winttcrlin, J., Trost, J., S, R., Schuster, R., Zambclli, T., Ertl, G.: Real- 
time STM observations of atomic equilibrium fluctuations in an adsorbate 
system: O/Ru(0001). Surf. Sci. 394 (1997) 159-169 

[45] Abramowitz, M., Stegun, I., eds.: Handbook of mathematical functions. 
National Bureau of Standards (1964) 

[46] Theil, F.: A proof of crystallization in two dimensions. Comm. Math. 
Phys. 262 (2006) 209-236 

[47] Aftalion, A., Serfaty, S.: Lowest Landau level approach in superconductiv- 
ity for the Abrikosov lattice close to H C2 . Selecta Math. 13 (2007) 183-202 

[48] Chen, X., Oshita, Y.: An application of the modular function in nonlocal 
variational problems. Arch. Ration. Mech. Anal. 186 (2007) 109-132 

[49] Massari, U.: Esistenza e regolarita dcllc ipcrsuperfice di curvatura media 
assegnata in R n . Arch. Rational Mech. Anal. 55 (1974) 357-382 

[50] Gilbarg, D., Trudinger, N.S.: Elliptic Partial Differential Equations of 
Second Order. Springer- Vcrlag, Berlin (1983) 

[51] Giusti, E.: Minimal surfaces and functions of bounded variation. Volume 80 
of Monographs in Mathematics. Birkhauser, Basel (1984) 

[52] Choksi, R., Sternberg, P.: On the first and second variations of a nonlocal 
isoperimetric problem. J. Rcinc Angew. Math. 611 (2007) 75-108 

[53] Kerner, B.S., Osipov, V.V.: Phenomena in active distributed systems. 
Mikroelektronika 14 (1985) 389-407 



47 



[54] Muratov, C.B.: Instabilities and disorder of the domain patterns in the 
systems with competing interactions. Phys. Rev. Lett. 78 (1997) 3149- 
3152 

[55] Muratov, C.B., Osipov, V.V.: General theory of instabilities for pattern 
with sharp interfaces in reaction-diffusion systems. Phys. Rev. E 53 (1996) 
3101-3116 

[56] Fusco, N., Maggi, F., Pratelli, A.: The sharp quantitative isoperimetric 
inequality Ann. of Math. 168 (2008) 941-980 

[57] Muratov, C.B.: Synchronization, chaos, and the breakdown of the collective 
domain oscillations in reaction-diffusion systems. Phys. Rev. E 55 (1997) 
1463-1477 

[58] Kiessling, M.K.H., Spohn, H.: A note on the eigenvalue density of random 
matrices. Comm. Math. Phys. 199 (1999) 683-695 

[59] Struwe, M.: Variational methods : applications to nonlinear partial differ- 
ential equations and Hamiltonian systems. Springer, Berlin (2000) 

[60] Kerner, B.S., Osipov, V.V.: Stochastically inhomogeneous structures in 
nonequilibrium systems. Sov. Phys. - JETP 52 (1980) 1122-1132 

[61] Mimura, M., Tabata, M., Hosono, Y.: Multiple solutions of two-point 
boundary value problems of Neumann type with a small parameter. SI AM 
J. Math. Anal. 11 (1980) 613-631 

[62] Modica, L., Mortola, S.: Un esempio di T^-convergenza. Boll. Un. Mat. 
Ital. B 14 (1977) 285-299 

[63] Modica, L.: The gradient theory of phase transitions and the minimal 
interface criterion. Arch. Rational Mech. Anal. 98 (1987) 123-142 

[64] Attouch, H., Buttazzo, G., Michaillc, G.: Variational analysis in Sobolcv 
and BV spaces. Society for Industrial and Applied Mathematics, Philadel- 
phia (2006) 

[65] Fife, P.C., McLeod, J.B.: The approach of solutions of nonlinear diffusion 
equations to traveling front solutions. Arch. Rat. Mech. Anal. 65 (1977) 
335-361 



48 



